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ABSTRACT 

We present the broadband X-ray properties of four of the most X-ray luminous (Lx > 10 45 erg s _1 in 
the 0.5-2 keV band) radio-quiet QSOs found in the ROSAT Bright Survey. This uniform sample class, 
which explores the extreme end of the QSO luminosity function, exhibits surprisingly homogenous X- 
ray spectral properties: a soft excess with an extremely smooth shape containing no obvious discrete 
features, a hard power law above 2 keV, and a weak narrow/barely resolved Fe Ka fluorescence line 
for the three high signal-to-noise ratio (S/N) spectra. The soft excess can be well fitted with only a 
soft power law. No signatures of warm or cold intrinsic absorbers are found. 

The Fe Ka centroids and the line widths indicate emission from neutral Fe (E — 6.4 keV) originating 
from cold material from distances of only a few light days or further out. The well-constrained 
equivalent widths (EW) of the neutral Fe lines are higher than expected from the X-ray Baldwin 
effect which has been only poorly constrained at very high luminosities. Taking into account our 
individual EW measurements, we show that the X-ray Baldwin effect flattens above Lx ~ 10 44 erg 
s _1 (2-10 keV band) where an almost constant (EW) of ^100 eV is found. 

We confirm the assumption of having very similar X-ray AGN properties when interpreting stacked 
X-ray spectra. Our stacked spectrum serves as a superb reference for the interpretation of low S /N 
spectra of radio-quiet QSOs with similar luminosities at higher redshifts routinely detected by XMM- 
Newton and Chandra surveys. 

Subject headings: galaxies: active - quasar: general - X-ray: individual (RBS 1055, RBS 320, RBS 
1124, RBS 1883) 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are classified by lumi- 
nosity and by phenomcnological criteria. The most lu- 
minous objects are the radio-loud quasars (quasi-stellar 
radio sources) and the radio-quiet QSOs (quasi-stellar 
objects). Radio-loud objects are characterized by rel- 
ativistic jets. Since the jet emission of the radio-loud 
objects is likely to dilute the X-ray spectrum from the 
central engine of the AGN, radio-quiet objects are more 
suitable for the study of X-ray emission associated with 
the accretion disk. 

To date, detailed information on the X-ray spectra of 
radio-quiet AGNs has largely been derived from observa- 
tions of nearby flux-bright objects. Spectroscopy of low 
luminosity (Lx < 10 44 erg s _1 ) radio-quiet Seyfert galax- 
ies has revealed a relatively complex multi-component 
X-ray spectrum. In the soft X-rays, a steep component 
(soft excess) has been observed in many cases. In the 
past this component has been interpreted as the tail of 
the thermal emission from the accretion disk or Comp- 
tonization of EUV accretion disk photons (Arnaud et 
al. 1985; Kawaguchi et al. 2001; Porquet et al. 2004). 
On the other hand, Crummy et al. (2006) and Ross & 
Fabian (2005) discuss ionized disk reflection as an expla- 
nation. The spectrum at soft X-rays can also be influ- 
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enced by a cold or warm absorber. Following this ob- 
servation, Gierlihski & Done (2004) consider absorption 
along the line of sight from relativistically outflowing, 
optically thin, warm gas as an explanation for the soft 
excess. At higher energies, X-ray spectra are dominated 
by a hard (r~2) power law, attributed to inverse Comp- 
ton emission from a hot, thin corona surrounding the 
accretion disk. Furthermore, Compton reflection of the 
primary flux by the disk or torus is inferred by the detec- 
tion of a 6.4 keV fluorescent Fe Ka line and the 7.1 keV 
Fe absorption edge in AGN spectra. The first narrow, 
unresolved Fe Ka line in an AGN was observed in OSO- 
8 spectra of Centaurus A by Mushotzky et al. (1978). 
Various XMM-Newton and Chandra observations have 
shown that a narrow (FWHM up to ~ 10 3 km s _1 due to 
Dopplcr broadening) Fe line component is virtually ubiq- 
uitous in Seyfert spectra and indicates emission from gas 
roughly commensurate with the optical Broad Line Re- 
gion (e.g., Nandra 2006; Shu et al. 2010) and/or emission 
from the putative Compton-thick, homogeneous molecu- 
lar torus invoked in Seyfert type I/II Unification schemes 
and thought to be located about light-years from the su- 
permassive black hole (e.g., Antonucci 1993). A relativis- 
tically broadened (FWHM~ 10 5 km s _1 ) component in- 
dicates emission from within a few i?sch (= 2GMbh/c 2 ) 
of the black hole where Mbh is the black hole mass (e.g., 
Tanaka et al. 1995; Fabian et al. 2002) but appears in 
only about half of Seyferts (Nandra et al. 2007) . 

There is evidence for the narrow Fe Ka equivalent 
width to be anti-correlated with X-ray luminosity, the 
so-called X-ray Baldwin effect. This effect has been doc- 
umented in several AGN samples, including Iwasawa & 
Taniguchi (1993) using Ginga, Nandra et al. (1997) us- 
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ing ASCA, and Page et al. (2004a, 2005) using XMM- 
Newton. However, the fitted correlation is mainly due to 
L2-10 < 10 45 erg s" 1 objects (Bianchi et al. 2007). 

Improving the signal-to- noise ratio (S/N) by stack- 
ing the X-ray spectra of high-redshift AGNs is a com- 
mon method for high-redshift AGNs. Streblyanska et 
al. (2005) stack the XMM-Newton spectra of AGNs in 
the Lockman Hole survey, which covers a large range 
of luminosities and rcdshifts. In the average spectra 
of both type I and type II objects they find significant 
Fe fluorescent emission, which seemed to be relativisti- 
cally blurred (type I AGNs: EW~600 eV, type II AGNs: 
EW~400 eV). In comparison, recent studies such as the 
XMM-Newton Wide Angle Survey by Corral et al. (2008) 
find only a narrow unresolved Fe line of EW~90 eV in 
the stacked spectrum of 606 type I AGNs. Mateos et 
al. (2010) determine the mean continuum shape of 487 
AGNs in the XMM-Newton Wide Angle Survey by stack- 
ing the spectra at different cosmic epochs. 

The i?0&4TBright Survey (RBS; Schwope et al. 2000) 
identifies the brightest (0.1-2.4 keV count rate >0.2ct 
s _1 , flux limit /o 5-2 ~ 10~ 12 erg cm -2 s -1 ) X-ray 
sources from the ROSAT All-Sky Survey (RASS) at high 
Galactic latitudes. X-ray luminous radio-quiet QSOs are 
rare locally and their spectral properties are so far poorly 
established. Due to their high X-ray flux, the RBS ob- 
jects form an ideal sample to obtain the most statistically 
significant X-ray spectra and allows us to evaluate di- 
rectly their X-ray spectral properties and the properties 
(a, EW) of their Fe Ka lines. In this paper, we use XMM- 
Newton to study the spectral properties of four objects 
that belong to the 12 most X-ray luminous radio-quiet 
QSOs found in the RBS with L .5-2 > 10 45 erg s" 1 . 
Such analysis of the most luminous RBS- QSOs will test 
the assumption of having very similar X-ray AGN prop- 
erties when stacking their X-ray spectra. The sample 
will serve as a medium- z reference for the interpretation 
of low S/N QSO spectra at high redshifts. In addition, 
the EW of the Fe Ka line will be directly accessible in 
the individual X-ray spectra. This will add several data 
points in the poorly studied high luminosity range of the 
X-ray Baldwin effect. 

The paper is organized as follows. In Section 2 we 
describe how we selected the sample. Section 3 gives 
black hole mass and Eddington ratio estimates from op- 
tical spectra. In Section 4 the X-ray data reduction is 
specified, and we present the analysis of the X-ray data 
for RBS 1055, RBS 320, RBS 1124, and RBS 1883 in 
Section 5. In Section 6 we derive the optical-to-X-ray 
spectral indices for our objects. Section 7 comments on 
the long-term flux variability of the RBS-QSOs. The 
construction and analysis of the averaged spectrum of 
the most X-ray luminous radio-quiet RBS-QSOs is given 
in Section 8. In Section 9 we discuss the results in the 
context of other studies, and we give our conclusions in 
Section 10. 

All quoted uncertainties represent 90% confidence in- 
tervals (Ax 2 = 2.706 for one parameter) unless otherwise 
stated. We adopt a standard flat ACDM cosmology with 
n M = 0.3, O a =0.7, H = 70kms- 1 Mpc- 1 (Spergel et 
al. 2003). 



2. CONSTRUCTION A SAMPLE OF THE MOST X-RAY 
LUMINOUS RADIO-QUIET RBS-QSOS 

Based on the RASS (Voges et al. 1999), Schwope et 
al. (2000) present a spectroscopic identification catalog 
of high-Galactic (\b\> 30°) latitude RASS sources with 
0.1-2.4 keV count rates above 0.2 s _1 . The identifica- 
tion program of 2042 RBS sources reaches a spectroscopic 
completeness ratio of more than 99.5%. 

We preselect objects that are classified as AGNs 
and sort them in descending order of their Galactic- 
absorption corrected 0.5-2 keV luminosity. BL Lac ob- 
jects are not considered for this study. Since we are inter- 
ested only in radio-quiet AGNs, we furthermore exclude 
radio-loud AGNs by their radio-to-optical flux density 
R > 10 (Kcllcrmann et al. 1989). The NASA/IPAC 
Extragalactic Database 4 (NED) is used to obtain radio 
detections and fluxes. In addition, we check the object's 
coordinates for entries in the radio catalogs NRAO VLA 
Sky Survey (NVSS; Condon et al. 1998), Sydney Univer- 
sity Molonglo Sky Survey (SUMSS; Mauch et al. 2003), 
and Faint Images of the Radio Sky at Twenty centimeters 
survey (FIRST, White et al. 1997). Whenever available 
we employ the given 5 GHz flux densities. However, in 
a few cases we have to extrapolate the 1.4 GHz (NVSS 
and FIRST) or the 843 MHz (SUMSS) flux density to 
the 5 GHz density by using the following relation: flux 
density oc v ai with a r = —0.5 (Kellermann et al. 1989). 
The optical flux densities are calculated from the U.S. 
Naval Observatory survey (USNO; Monet et al. 1998) 
.B-band magnitude at an effective wavelength of ~4100 
A (Gajdoski & Weinberger 1997). 

Introducing a low luminosity cut of £0.5-2 = 10 45 erg 
s , we construct a sample of the 14 most X-ray luminous 
radio-quiet RBS-QSOs. Table 1 shows the main proper- 
ties of the sample. The objects are ordered in descend- 
ing X-ray luminosity and count rate in the 0.5-2 keV 
ROSAT band. We list the position of the optical coun- 
terpart (taken from Schwope et al. 2000), the USNO B- 
and i?-band magnitude, the RASS 0.5-2 keV luminosity 
(corrected for Galactic absorption), the RASS count rate 
in the hard R OS A T band (0.5-2 keV), the redshift of the 
optical counterpart, the 5 GHz to B-band flux density 
ratio R (a value of R = means that the object was not 
detected as a radio source) , and the instrument which ob- 
served the corresponding object. All RBS objects listed 
in Table 1 have optical counterparts with broad emission 
lines in their spectra (type I AGN) . The redshifts of our 
RBS sources are based on the values listed in the NED, 
since it contains the latest available information on the 
objects. 

The first two most X-ray luminous radio-quiet RBS- 
QSOs are RBS 1423 with L Q 5 _ 2 = 10 46 5 erg s" 1 and 
RBS 825 (Lo.5-2 = 10 46 1 erg s" 1 ). However, RBS 1423, 
whose identification with a z — 2.26 QSO was somewhat 
uncertain, indeed was identified using the XMM-Newton 
position with a low luminosity (£0.5-2 = 2.1 x 10 44 erg 
s _1 ) z = 0.208 QSO with a significantly broadened ion- 
ized Fe Ka line (Krumpe et al. 2007). Object RBS 825 
turned out to be a quadruple QSO lens system with 
one of the largest known separations (Inada et al. 2003). 
Hence, its observed flux is significantly enhanced by lens- 
ing and its intrinsic luminosity much lower (Lamer et al. 
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Table 1. Observed Properties of the Most X-ray Luminous Radio-quiet 
QSOs in the ROSAT Bright Survey 



Object R.A. Decl. B band R band log^x/Iergs -1 ]) RASS CR z R Observations 



Name 


(J2000) 


(J2000) 


(mag) 


(mag) 


(0.5-2 kcV) 


(0.5-2 keV) F 5 GHjF 0pt 




RBS 1423 3 - 


14 44 14.7 


+06 33 07 


16.3 


16.7 


46.5 


0.11 


0.208 





XMM 


RBS 825 h 


10 04 34.9 


+41 12 40 


18.0 


18.2 


46.1 


0.08 


1.73 





XMM 


RBS 1055 


11 59 41.0 


-19 59 25 


16.5 


16.1 


45.3 


0.19 


0.450 


5 


XMM 


RBS 1273 


13 29 28.6 


-05 31 36 


15.3 


15.4 


45.2 


0.09 


0.580 







RBS 320 


02 28 15.2 


-40 57 11 


15.4 


14.9 


45.1 


0.10 


0.494 





XMM 


RBS 1124 


12 31 36.5 


+70 44 14 


15.3 


15.4 


45.0 


0.45 


0.209 





XMM/Suzaku 


RBS 229 


01 40 17.0 


-00 50 03 


16.5 


16.3 


45.0 


0.17 


0.334 


1 




RBS 832 


10 08 31.5 


+46 29 53 


19.3 


18.5 


45.0 


0.14 


0.388 







RBS 1062 


12 03 43.2 


+28 35 57 


17.2 


17.7 


45.0 


0.12 


0.374 







RBS 1502 


15 28 40.7 


+28 25 29 


16.1 


16.6 


45.0 


0.08 


0.450 







RBS 1853 


22 20 49.7 


-31 56 54 


17.6 


17.2 


45.0 


0.08 


0.503 







RBS 1883 


22 41 55.3 


-44 04 58 


15.1 


15.7 


45.0 


0.06 


0.545 





XMM 


RBS 2040 


23 43 13.5 


-36 37 54 


15.4 


15.8 


45.0 


0.05 


0.620 


2 




RBS 218 


01 35 20.2 


-29 23 58 


16.7 


15.9 


45.0 


0.04 


0.699 








Notes. The four objects in bold face are the sources discussed in this paper. 
a Corrected; no longer among the most X-ray luminous (AGN at z = 0.208 instead of an AGN at z = 2.262, log(Lx/[crg s -1 ]) = 44.3). 
b Corrected; no longer among the most X-ray luminous (lensed system). 



Table 2. Fit Results from Optical Spectra, Black Hole Masses, and 
Eddington Ratios for the Observed RBS-QSOs 



Object 
Name 


Line 


FWHM 

(km s- 1 ) 


A 

(rcst-f.) (A) 


(erg s" 1 ) 


Reference 


log(M B H/M Q ) 


(10 44 erg s- 1 ) 


ibol/^edd 


RBS 1055 


Mg II 


5540 


3000 


2.00X10 45 


Woo (2008) 


8.87 


118 


0.13 


RBS 320 


H/9 


2265+100 


5100 


1.21 xlO 46 


Grupc ct al. (1999) 


8.88 


1370 


1.43 


RBS 1124 


H/9 


4260+1250 


5100 


1.26X10 44 


Grupe et al. (2004) 


8.06 


14.2 


0.10 


RBS 1883 


H/9 


1890+200 


5100 


2.58X10 45 


Grupc ct al. (1999) 


8.26 


292 


1.27 



2006). 

Removing the exceptional objects RBS 1423 and RBS 
825, the most X-ray luminous radio-quiet RBS-QSO RBS 
1055 is only the 54th most X-ray luminous object within 
all RBS entries (20th most X-ray luminous when exclud- 
ing BL Lac objects). Rephrasing this interesting result, 
no radio-quiet QSO with log(Lo.5-2/[ergs -1 ]) > 45.3 is 
found in this flux limited all-sky survey. In contrast, the 
three most X-ray luminous radio-loud RBS-QSOs (RBS 
717, RBS 315, and RBS 1788) have L .5-2 = 10 4> erg 
s _1 and R > 2000. This suggests that the jet in radio- 
loud QSOs or jet-related changes in the accretion disk 
significantly boost the X-ray luminosity. 

We have observed two of the most X-ray luminous 
radio-quiet RBS-QSOs which have high 0.5-2 keV RASS 
count rates (RBS 1055 and RBS 1124). We combine 
those with two archival data sets (RBS 320 and RBS 
1883) to create a subsample of four objects that are dis- 
cussed in this paper. 

3. BLACK HOLE MASS ESTIMATES FROM OPTICAL 
SPECTRA 

All four RBS-QSOs analyzed in this study have been 
spectroscopically observed in the optical. Over the last 
decade reverberation mapping of AGN broad emission 
lines established a simple empirical method to estimate 
black hole masses of a large population of AGN by mea- 
suring the FWHM and the luminosity which are acces- 
sible in a single spectrum (size-luminosity relation; see, 
e.g., Kaspi et al. 2000; Vestergaard & Peterson 2006). In 
this section, we give mass estimates of the central black 



hole and the Eddington ratio -Lboi/^cdd for each of our 
objects based on published reverberation mapping (RBS 
1055, also known as CTS 306; Woo 2008) and published 
line property and luminosity measurements from optical 
spectra (RBS 320, RBS 1124, RBS 1883). 

Table 2 lists the FWHM values and the related contin- 
uum luminosity which are used to infer black hole masses, 
the reference of these values, and the derived quantities. 
Grupe et al. (1999) give only the values for FWHM of H/3, 
but not the XL\ continuum luminosity. We used their op- 
tical spectra of RBS 320 and RBS 1883 (Beuermann et al. 
1999; H.-C. Thomas, private communication) to compute 
the relevant fluxes which were corrected for Galactic ex- 
tinction. The Mbh estimates are based on the empirical 
formula from McGill et al. (2008). The bolometric lumi- 
nosities are estimated by multiplying a factor of 5.9 and 
11.3 (McLure & Dunlop 2004) to the 3000 A and 5100 A 
continuum luminosity, respectively. The Eddington ra- 
tio Lhoi/^edd is the ratio between the bolometric and the 
Eddington luminosity (L odd = 1.26 x 1O 38 M B h/M erg 
s^ 1 ). Note that Grupe et al. (2004) estimate the bolo- 
metric luminosity from a combined fit to optical-UV and 
X-ray data to be Lboi = 3.4 x 10 45 erg s _1 . Grupe et al. 
(2010) list updated values for RBS 1124 based on simul- 
taneous optical-to-X-ray observations with Swift and find 
an accretion rate relative to Eddington Lb o i/£cdd=0.63. 
However, we use the correction factor to derive an esti- 
mate of the bolometric luminosity to treat all sources in 
a consistent way. 

The values for Mbh, -^boi, and Lboi/^edd in Table 2 are 
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Table 3. XMM-Newton Observations of the Most Luminous Radio-quiet 
QSOs in the ROSAT Bright Survey 



Object 
Name 


Date 
(Year) 
(MM-DD) 


ObsID 


Nh,gal 
(10 20 cm" 2 ) 


Wl\/T 
VV 1V1 

pn 
MOS 


Filter 

pn 
MOS 


IN CI ri*xp . 

pn (ks) 
MOS (ks) 


fP .r noioii 
pn 
MOS 


Count Rate 
pn (ct s -1 ) 
MOS (ct s- 1 ) 


/X 

(0.5-2 keV) 
(2-12 keV) 


(0.5-2 keV) 
(2-12 keV) 


RBS 1055 


2008 


0555020201 


4.12 


PFW 


T 


21.3 


42320 


1.94±0.01 


2.0 


15 




06-11 






PFW 


T 


24.8 


13904 


0.545±0.005 


4.1 


28 


RBS 320 


2004 


0200480101 


2.20 


PFW 


T 


26.7 


42533 


1.57±0.01 


1.2 


16 




12-07 






PFW 


T 


31.3 


10833 


0.339±0.003 


0.93 


9.6 


RBS 1124 


2008 


0555020101 


1.69 


PFW 


T 


15.5 


54901 


5.36±0.02 


4.6 


5.8 




05-15 






PPW3 


T 


21.1 


28701 


1.44±0.01 


7.1 


8.9 


RBS 1883 


2003 


0153220101 


1.82 


PFW 


M 


0.5 


422 


0.70±0.04 


0.57 


11 




05-17 






PFW 


H 


5.2 


670 


0.118±0.005 


0.40 


5.6 



Abbreviations: WM - window mode of the different detectors, PFW - prime full window mode, PPW3 - prime partial W3 mode (to 
decrease the risk of photon pile-up), T - thin filter, M - medium filter, H - thick filter, Net Exp. - net exposure time (average MOS 
time for each MOS detector), # Photon - source photon number in the 0.2-12 keV band. Count rates in the energy band of 0.2-12 keV. 
Intrinsic, Galactic absorption-corrected fluxes (in units of 10~ 12 erg enr 2 s _1 ) and rest-frame luminosities (in units of 10 44 erg s _1 ) arc 
based on the best-fit model of each object. 



subject to large uncertainties and should be understood 
as approximate quantities only. The super-Eddington 
ratios of RBS 320 and RBS 1883 prompt one to investi- 
gate the properties of the H/3 lines in more detail. RBS 
320 has to be modeled with a broad and narrow Gaus- 
sian line profile. Considering only the broad line profile 
with an instrumentally corrected FWHM= 2900 ± 200 
km s^ 1 leads to an upper limit of the black hole mass 
log(MBH/Af©) < 9.09 and a lower estimate of the Ed- 
dington ratio Lboi/^cdd >0.88. In the case of RBS 1883, 
the red wing of H/3 falls in the atmospheric A-band mak- 
ing the modeling of the line even more difficult. Cor- 
recting for the ^4-band absorption, we only needed one 
broad line Gaussian profile with an instrumentally cor- 
rected FWHM= 3000 ± 300 km s" 1 to fit the shape of 
the H/3 line. This results in log(M B H/M Q ) = 8.66 and 
an Eddington ratio £hoi/£cdd=0.51. 

Despite all systematics and uncertainties in the esti- 
mate, the extremely high X-ray luminosities in these 
radio-quiet QSOs originate in a combination of very high 
black hole masses and very high Eddington ratio, but 
they are not extreme (e.g., Schulze & Wisotzki 2010; 
Steinhardt & Elvis 2010). 

4. X-RAY DATA REDUCTION 

All XMM-Newton data are processed with SAS version 
7.1.0 (Science Analysis Software, Gabriel et al. 2004) 
package, including the corresponding calibration files. 
The tasks emproc and epproc are used for generating 
linearized photon event list from the raw EPIC data. 
The RGS data are reduced by the task rgsproc. The 
data are cleaned for times of high particle backgrounds. 
We investigate the astrometrically corrected EPIC coor- 
dinates to verify the corresponding optical counterpart. 
All identifications are confirmed by our XMM-Newton 
data. 

We make use of the Pipeline Processing System (PPS) 
light curves. The light curves are background, dead time, 
and mirror vignetting corrected and are extracted in an 
energy range of 0.2-12 keV. We rebin the light curves 
to bins of 1000 s. For RBS 1883 only the MOS light 
curves are considered, because of the high background 
count rate in the pn detector. 

The X-ray spectra are extracted with the task 



especget. For the pn detector we use box-shaped back- 
ground regions at roughly the same detector y-position 
as the circular source region is extracted. The pn back- 
ground regions are located on the same chip as the source 
region, except for RBS 1124 where we are forced to use 
the adjacent chip due to the large covered area of RBS 
1124's source region. For the MOS detectors we use 
annulus-shaped background regions around the source 
on the same CCD. Additional sources in the background 
regions are excluded. Furthermore, we follow the rec- 
ommended flag selection of the macros XMMEA_EP and 
XMMEA_EM. The low energy cut-off is set to 0.2 keV, and 
the high energy cut off is fixed to 12 keV. Using the task 
epatplot, we verify that the effect of photon pile- up is 
negligible in all sources, except RBS 1124. We select 
X-ray events corresponding to pattern 0-12 (single and 
doubles) for the MOS detectors and 0-4 events for the 
pn detector, except for RBS 1124 where we only select 
single events for the pn detector. The count rates of 
RBS 1124 (Table 3), which is the flux-brightest radio- 
quiet RBS-QSO, are close to the count rate where pho- 
ton pile- up is expected (pn: ^6 ct s _1 ; MOS operated 
in large window-prime partial w3 mode: ~1.8 ct s _1 ; 
XMM-Newton Users Handbook 5 ). The MOS detectors 
are not affected by photon pile-up. However, the pn de- 
tector shows signs of pile-up. Therefore, we follow the 
recommended procedure of defining the source region as 
an annulus to exclude the central region (r = 2'.'5) with 
the highest pile-up probability. Furthermore, we only 
use the pn pattern spectrum which is less sensitive to 
pile-up than the pattern 1-4 spectrum. The resulting 
0.2-12 keV count rate for the pn detector is therefore 



cr pn = 3.49 ±0.02 ct s" 1 . 



The spectra are binned to a minimum of 20 counts per 
bin to apply the ^-minimization technique. The result- 
ing X-ray spectra are fitted with Xspec version 12.5.0 
(Arnaud 1996). We use the cosmic abundances accord- 
ing to Wilms, Allen, & McCray (2000) and the photo- 
electric absorption cross sections provided by Verner ct 
al. (1996). The Galactic hydrogen column densities in 
the line of sight to the objects (Table 3) are derived from 

5 http : // xmm . esac . esa. int/external/xmm_user_support/documentat ion/uhb 
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the -/Vn-calculator 6 which is based on Dickey & Lockman 
(1990). 

4.1. RBS-QSOs Light Curves 

To investigate possible variability of the sources which 
may be related to spectral variability, we first study the 
PPS light curves of our RBS-QSO sources. The visual 
inspections verify that for RBS 1055, RBS 1124, and RBS 
1883 no significant variability is present. 

In the case of RBS 320 the average count rate stays 
nearly constant up to the middle of the observation and 
then increases almost monotonically. At the end of the 
observation the count rate is ^20% higher than the av- 
erage count rate at the beginning of the observation (the 
background is not showing this trend). The statistical 
ltr fluctuations are on the order of 3%. 

To investigate spectral variability, we extract two spec- 
tra for RBS 320, one covering the first two-thirds of 
the observation, the other one the last third. We ap- 
ply the best-fit time-average model to each sub-spectrum 
(2 power laws + Gaussian line, 2PLG; see Section 5.2), 
keeping as many parameters frozen at the time-average 
values as possible. We first thaw only the normalization 
of the hard power law. x 2 /d.o.f. is 714/693 and 522/501 
for the low and high flux spectra, respectively. When 
we thaw the normalization of the soft power law as well, 
X 2 /d.o.f. falls in each case to 710/692 and 517/500. Con- 
sequently, the change of the normalization for the hard 
power law is more significant than for the soft power law 
when considering the uncertainties. The Fe line parame- 
ters are consistent with their time averaged values. Fur- 
thermore, there is no significant evidence for variations 
in the photon index of either power law. We conclude 
that the increase in the count rate of RBS 320 toward 
the end of the observation is consistent with a change in 
power-law normalization mainly in the hard component, 
consistent with the notion that the physical properties 
of the corona (e.g., temperature) are likely not chang- 
ing. The analysis of the time-averaged spectrum over 
the whole observation is therefore justified. 

5. SPECTRAL ANALYSIS 

5.1. X-ray Spectra of RBS 1055 

We start with the simplest modeling of the 0.2-12 
keV data of the most X-ray luminous radio-quiet RBS- 
QSO: a simple power law (henceforth called the "1PL" 
model); this and all models discussed herein are corrected 
for Galactic absorption (see Table 3). For all objects, 
studied in this paper, no significant evidence for neu- 
tral absorption in excess of the Galactic column den- 
sity is found. As the 1PL model is a poor fit (x 2 /d.o.f. 
= 1413/1251), we then try two power laws, one mod- 
eling the soft, the other the hard energy range (hence- 
forth called the "2PL" model). x 2 /d.o.f. improves to 
1288/1249. The only residuals left are positive residuals 
around 6.4 keV in rest frame (4.4 keV in observed frame 
as shown in Figure 1). We cannot successfully model 
the soft excess using a blackbody component (hence- 
forth called the "PLBB" model; x 2 /d.o.f. = 1319/1249, 
k B T = 0.11 ± 0.01 keV, V = 1.73 ± 0.01). We return 
to the 2PL and add a Gaussian fixed at 6.4 keV (rest 

6 heasarc . gsf c .nasa. gov/cgi-bin/Tools/w3nh/w3nh .pi 




observed energy (keV) 

Fig. 1.— XMM-Newton EPIC spectra of RBS 1055 in the ob- 
served frame. The spectra are fitted with a single power law (1PL) 
to data in the 1—12 keV range. At energies below 0.5 keV a soft 
excess is visible. A fluorescent 6.4 keV Ee Ka line will be shifted 
to 4.4 keV at z = 0.450. 



frame) consistent with neutral Fe Ka emission (hence- 
forth called the "2PLG" model). The fit to the data 
(x 2 /d.o.f. = 1266/1247) and the residuals are signifi- 
cantly improved. 

The shape of the spectrum near the Fe line cannot 
be explained with alternative models incorporating ab- 
sorption. For example, we try a partial covering model 
incorporating absorption by gas with a column density of 
Ah ~ 6 x 10 23 cm~ 2 and a covering fraction of ^30%, but 
this fails to properly model the residuals near 3-6 keV 
(x 2 /d.o.f. = 1274/1249). Although there are no obvi- 
ous strong absorption lines or edges (e.g., O VII, O VIII) 
present in the spectrum, we also test full-covering and 
partial-covering absorption by ionized material. We use 
an XSTAR 7 table (turbulent velocity v = 200 km s -1 , in- 
put continuum with a photon index r = 2), but in either 
case, there is no improvement to the fit and the data 
points in the 3-6 keV band are not properly modeled. 
For the full-covering absorption model by ionized mate- 
rial the 2rj upper limit estimate for the column density is 
A H , warm = 3.9 x 10 21 cm" 2 for log (f /(erg cm s- 1 ))^ 1.7. 
In the case of the partial-covering model the 2a upper 
limit is Ah, warm = 1-6 x 10 21 cm~ 2 for log (£/(erg cm 
s _1 ))= 1.7. In addition, no strong Fe K edge at 7.1 
keV is found. Adding a zedge component to the 2PL 
model with an energy fixed at 7.1 keV, yields an upper 
limit of t — 0.3 and does not improve the fit (x 2 /d.o.f. 
= 1288/1248). All these models do not account for the 
positive residuals around 6.4 keV as well as the 2PLG 
model does. Consequently, we conclude that the Fe Ka 
emission line is real. 

Leaving the line energy as a free parameter slightly im- 
proves the fit (x 2 /d.o.f. = 1264/1246; henceforth called 
the "2PLG2" model). We find E FCtIcst = 6.25 ± 0.22 

keV and <7Fc,rcst = 0.49to'is keV. Using the response 
matrix, we compute an instrumental energy resolution 
of FWHM re soi ~ 160 eV at the line energy. For all ob- 
jects we will only consider a line to be significantly re- 
solved if the line is resolved at the 95% (2a) confidence 
level. Since a fit with an unresolved component (Gaus- 
sian line profile with a line centroid fixed at 6.4 keV and 
width a = 0.001 keV) results in a worse fit (x 2 /d.o.f. = 

7 http : //heasarc. gsf c .nasa. gov/docs/sof tware/xstar/xstar . html 
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1277/1247) and the lower limit 95% confidence errors of 
the measured <XFe,rest does not include a = 0.001 keV, 
we conclude that we have detected a resolved, neutral 
Fe Ka line with EW= 0.14±°;^ keV (3er uncertainties) 
in the spectrum of RBS 1055. We exclude the presence 
of an additional, unresolved component: we add a sec- 
ond narrow Gaussian with energy centroid fixed at 6.4 
keV (rest frame) and a fixed at 1 cV, but this yields no 
improvement to the fit. 

The freely fitted line energy of E = 6.25 keV and the 
broad shape could potentially be caused by the presence 
of a Compton shoulder to the 6.4 keV Fe line. If so, we 
expect a Compton shoulder intensity of 20%-40% of the 
central Fe line, with a centroid energy of ^6.3 keV (Matt 
2002). Fixing the line energy of the primary Gaussian 
profile to 6.4 keV and including an additional Gaussian 
profile at a fixed energy of E = 6.3 keV with an intensity 
of 30% of the primary Gaussian and a a = 0.07 keV 
(Bianchi ct al. 2002) does not improve the fit. Leaving 
the line intensity as a free parameter causes its intensity 
to tend to zero and making the primary Gaussian the 
only significant source for modeling the Fe line. 

We then tried to model the observed Fe line profile 
using a laor component (accretion disk around a Kerr 
black hole; Laor 1991), with an emissivity index fixed to 
3, inclination to 30°, and the outer radius to 400 gravita- 
tional radii (r g = GMbh/c 2 ). We achieved a fit virtually 
identical to that using a single Gaussian emission line 
(see above). The inner radius of emission is I6I7 6 r g> i- e > 
poorly constrained due to the weakness of the line. Sim- 
ilar results are found for the high S/N spectra for RBS 
320 and RBS 1124, and this model is thus not discussed 
further. 

For RBS 1055, 320, and 1124, we determined the 90% 
confidence level upper limits of the equivalent width of 
a relativistically broadened Fe line in addition to a nar- 
row profile by first modeling the Fe line with a Gaussian 
component (width a fixed to 1 eV) and then adding a 
laor component, leaving only the normalization of the 
laor component as a free parameter (R- m = 1.235 r g , 
Ri n = 400 r g , and emissivity index 3). The resulting 
upper limits are discussed in Section 9.3. 

To test if the soft excess and the broad Fe Ka line can 
be modeled in a self-consistent manner, we employed the 
ionized reflection model ref lion. We use a model of 
the form (kdblur*ref lion)+zpo (Ross & Fabian 2005). 
Given the limited S/N of our spectra, some model pa- 
rameters have to be fixed: we constrain the inclination 
to be 30°, the outer radius of the reflecting disk to 400 
gravitational radii, the emissivity index q to 3, and the 
photon index of the reflected power law equal to that of 
the illuminating power law. Two almost identical x 2 ~fit 
minima are found (the reduced \ 2 values are very sim- 
ilar): one fit (REFL1) is driven by both the soft excess 
and the Fe line and results in a low value of the ion- 
ization parameter; the second (REFL2) is driven by the 
shape of the soft excess and results in a high value of the 
ionization parameter. 

For REFL1, the best-fit (x 2 /d.o.f. = 1272/1247) 
yields an Fe abundance consistent with solar (Fe/solar 
= 1.4+Q g), an inner radius of R m = 16. OT^^ r g , and 
the lowest possible value of the ionization parameter, 
Cdisk = 30t „ erg cm s _1 (best-fit value pegged at hard 



limit of model); T = 1.77 ± 0.01. The fraction of the 
reflected component to the total 0.2-12 keV output is 
9%. The low ionization parameter models the Fe line 
well, but also causes relatively sharp line features in the 
soft excess model. Although a higher intensity of the re- 
flected component would still be in agreement with the 
modeling of the observed Fe line, it would also increase 
the sharp line features in the soft excess which are not 
present in the data. Therefore, the best-fit is an attempt 
to explain both features simultaneously. 

In the best-fit to REFL2, the fit (x 2 /d.o.f. = 
1268/1247) finds Fe/solar = 0.2^%^, an inner radius 
of R in < 6.3 r g , and £ disk = 15001^° erg cm s _1 ; 
r = 1.601 

o!o9- The reflection model comprises 34% of the 
total 0.2-12 keV output and the parameters are driven to 
describe the smooth soft excess in the data up to 3 keV. 
This reflection component has only a very minor contri- 
bution to the high energies. Due to the high ionization 
parameter, the Fe Ka line contribution is very broad, 
smeared out, and the line intensity too low. A higher Fe 
abundance is needed to improve the modeling of the Fe 
K line profile. However, this would cause a change in the 
shape of the soft excess around 0.7-0.9 keV in the rest 
frame, the energy range associated with Fe L emission, 
because the Fe K and Fe L line intensities are connected. 
Consequently, this model is also a compromise between 
accounting for the soft excess and representing the Fe 
line. Ignoring the energies around the Fe line still results 
in two equally good fit minima except that both are now 
found with the lowest Fe abundances possible. 

Leaving the emissivity index q free to vary in both 
models leads to a somewhat unconstrained emissivity in- 
dex and improves the fit only marginally. We keep the 
values of the emissivity index frozen at 3 unless a free in- 
dex results in improving the fit significantly (see Table 5). 
These values may not necessarily be representative of the 
physical/intrinsic emissivity indices, but such values are 
required to yield a soft excess with a smooth shape and 
thus achieve a good fit to the data. 

The blurred reflection model can only successfully 
model Fe lines that originate in close proximity to the su- 
permassive black hole and are relativistically broadened 
features. Therefore, we consider an additional narrow Fe 
line component with E fixed at 6.4 keV and a fixed at 
0.001 keV to model emission from gas located at large 
distances from the central region. This yields REFL1G 
and REFL2G, which differ from REFL1 and REFL2, re- 
spectively, by the addition of the narrow Gaussian line 
component. The residuals around 6.4 keV (rest frame) 
are better modeled in both cases (up to A\ 2 = —6; Ta- 
ble 5), and we consider REFL2G to be the best-fit model 
for this spectrum. 

We also try to incorporate two reflection components 
to account for reflection in different region: one very close 
to the central engine producing the smooth shape of the 
soft excess, the other far away to produce the narrow Fe 
line. However, the moderate quality of the XMM-Newton 
data for all objects does not allow us to constrain crucial 
parameters such as the break radius (radius of inner and 
outer reflection) and ionization parameter of the inner 
disk. 

The RGS data of RBS 1055, as well as for all other 
high S/N RBS-QSOs, are well fitted with a single power 
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Table 4. Spectral Power-Law Fits for RBS 1055, RBS 320, RBS 1124, 
and RBS 1883 



Model 


r 


^Fe.rest 

(keV) 


Te.rest 

(eV) 


EWpe.rest 

(eV) 


X 2 /d.o.i. 


IPL 
2PL 
2PLG 
2PLG2 


1.79±0.01 
1.3718$ & 2.121°;^ 
1.131°:?? & 1.98l8;? 

i.o7i8;2| & i.97±g;SS 


JrLh5o 1U55 

!6.4! 
6.25 ±0.22 


4801?^ 
490l?« 


14011° 
140l«8 


1413/1251 
1281/1249 
1266/1247 
1264/1246 


IPL 
2PL 
2PLG 
2PLG2 


2.69±0.01 

i.9ii8;» & 3.3118;!? 

1.95±0.10 & 3.3218;^ 
1.95l8;lJ & 3.32±0.12 


!6.4! 

6 47+ - 21 
D -^'-0.22 


270+ 420 
260+ 370 


120+ 7 J° 
13011° 


1615/810 
934/808 
924/806 
924/805 


IPL 
2PL 
2PLG 
2PLG2 


1.96±0.01 
1.13 ±0.04 & 2.19 ±0.03 
1.15 ± 0.05 & 2.19 ± 0.04 
1.15 ±0.05 & 2.19 ±0.04 


RBS 1124 

!6.4! 
6.49 ±0.14 


~~-(-320 

180l? 2 8 

140+240 
lw -140 


_,„+60 
70 

,u -40 
7Q +60 
' u -40 


1660/1293 
1525/1291 
1516/1289 
1515/1288 


IPL 
2PL 
2PLG 
2PLG2 


1.79±0.01 

i ro+0.08 & g ~,+0.29 
1 - os -0.10 '4-0.30 

1.58 ±0.08 & 2.68t° f 3 
i.eotlil &2.74±°f s 


!6.40! 

c , n +0.0S 


/ 20+ 140 


r ? +34 
01 -26 
rp+34 


1122/926 
958/924 
942/922 
942/921 


IPL 
2PL 


2.9±0.1 
1.7±0.3 & 3.9l8 j 


RBS 1883 






180/86 
119/84 



Note. — Models (explanation of models according to the Xspec notation): IPL - single power law (phabs (zpo) ), 
2PL - two power laws, one fitting the soft energy range, the other, the hard energy range (phabs (zpol+zpo2)), 
2PLG - two power laws including a Gaussian component with a line energy fixed to 6.4 keV in the rest frame 
(phabs (zpol+zpo2+zgauss) ), 2PLG2 - same as 2PLG but line energy is allowed to vary. Frozen parameters are in 
exclamation marks. Energies and equivalent widths are given in the quasar rest frame. Uncertainties are computed 
for the 90% confidence interval. For RBS 1124 we show the results of the 0.35-10 keV Suzaku data fits in italic font, 
while all 0.2-12 keV XMM-Newton-hascA fits are in roman font. 



Table 5. Results for REFL1G and REFL2G: Disk Reflection Plus a 
Narrow Fe Line component 



X 2 /d.o.f. 


r 


Q 


Rin 

(rg) 


(erg cm s x ) 


Fe/Solar 


Ax 2 
Fe Line 


1271/1246 
1262/1246 


-1 77 +0.02 
' -0.01 

1 60+ 07 
1 - ou -0.09 


!3! 
!3! 


RBS 1055 

15-911^ 
2-6l?L 


3013. 
1600l 7 88 


1 4+0-5 

o.2l8:?. 


1 

6 


933/804 
971/804 


2 fi[ -+0.08 

9 oq+0.02 
z - oa -0.03 


8.4l 2 ;f 

7.2l : 9 


RBS 320 

r, O + 0.3 

z -°-0.4 
1 O+0.6 
1 - a -0.3 


no±S82 

1500 ± 200 


0.1I8; 2 
0.3 ±0.1 


7 
10 


1541/1288 
1514/1288 


1-3 -0.02 
1 n r+0.09 


!3! 
!3! 


RBS 1124 

1+2.5 

■ i - 1 -1.9* 

7 r, + 2.8* 

'• z -4.4 


3018. 
9700+ 300 * 

yluu _1600 


3-5l 2 2 i 
2.8l?;8 


1 

2 


930/921 
931/920 


1 17 +0.10 
1 ' 1 '-0.06 


!3! 

s.otli* 


J -*-o* 

y 9+0.4 


, n +S4 

7000 +3000* 
iuuu_ 500Q 


, 7+0.9 
1 - '-0.8 

3.3±H 


4 

14 



Note. — For explanation see Table 4. A lower/upper 90% confidence uncertainty marked with "*" symbolizes that 
the parameter pegged at the lower/upper hard limit of the model. All fits use frozen parameters of: inclination angle 
= 30°, .Rout = 400 r g , £ F c,rcst = 6.4 keV, and er Fo , rcst = 1 eV. 



law (x 2 /d.o.f. = 207/240). No significant absorption 
or emission features are visible based on only the RGS 
data of RBS 1055. The upper limits to the absolute val- 
ues of equivalent widths (90% confidence level) of typical 
prominent narrow emission and absorption lines such as 



O VII, O VIII, Ne IX and Ne X are as high as ~ 8 eV. 
Both disk reflection scenarios fit equally well the RGS 
data. 

Considering the best-fit model for RBS 1055 
(REFL2G), we give the intrinsic (corrected for Galac- 
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tic absorption) 0.5-2 keV and the 2-10 keV fluxes of the 
object as well as the corresponding intrinsic luminosities 
for all objects in Table 3. 

5.2. X-ray Spectra of RBS 320 

For the analysis of the RBS 320 data (and the re- 
maining RBS objects), we follow the same approach as 
described in detail for RBS 1055. We will only sum- 
marize the most important findings. As for RBS 1055, 
we also find that a 2PL model is required to character- 
ize the overall shape of the RBS 320 data appropriately 
(x 2 /d.o.f. = 934/808) although positive residuals are 
still present at 6-7 keV (rest frame) where the Fe Ka 
line is expected (see Figure 2). 

A PLBB model (k B T = 0.09 ± 0.01 keV, T = 2.42 ± 
0.02) cannot fit the data well (x 2 /d.o.f. = 1099/808). 
The best-fit partial covering model incorporating absorp- 
tion by a gas with a column density of Ah ~ 2 x 10 
cm~ 2 and a covering fraction of ~60% fails to properly 
model the data (x 2 /d.o.f. = 1071/808). A moderately- 
good fit is obtained for full-covering absorption model 
with ionized material (x 2 /d.o.f. = 935/806). However, 
this fit finds a very high ionization value (log (£/(erg cm 
s _1 ))= 6, lower 90% limit: log (£/(erg cm s" 1 ))^ 3.3) 
and an unconstrained Ah, warm value. Ionization param- 
eters log (£/(erg cm s _1 ))> 4 do not change the overall 
shape of the X-ray spectrum; they only include narrow 
H-like absorption lines which are not seen in the data. 
Therefore, we also exclude for RBS 320 the presence of 
a warm absorber. An edge component with energy fixed 
at 7.1 keV does not improve the fit of the 2PL model 
(x 2 /d.o.f. = 934/807, 90% confidence upper limit on 
t = 0.19). Consequently, also in the case of RBS 320 we 
conclude that an Fe Ka emission line is detected. 

Adding a Gaussian component to the 2PL model re- 
sults in an improved fit (2PLG, x 2 /d.o.f. = 924/806, 
Table 4) and yields er F c,rcst = 0.271^ keV for an Fe Ka 
line with energy centroid fixed at E = 6.4 keV. The line 
is detected at more than 3cr confidence. Fixing a = 0.001 
keV results in a fit with A% 2 = +3. Following our general 
guideline for the detection of a resolved line, the line is 
not resolved. Leaving the line energy as a free parameter 
(2PLG2) does not change the fit. A Compton shoulder 
to the Fe line is not consistent with the data. 

Applying REFL1 with an emissivity index fixed to 3 
results in a moderately good representation of the data 
(x 2 /d.o.f. = 1032/806) when the lowest possible ioniza- 
tion parameter is used (£ = 30 to* erg cm s _1 ). The fit 
significantly improves when thawing the emissivity index 
(q = %tl*, X 2 /d.o.f. = 940/805). At q = 3 the soft ex- 
cess is still modeled with some structure not seen in the 
data and peaks at around 0.6 keV (rest frame; O VII), 
while much higher q values again result in a much more 
smooth soft excess without any structure and shift its 
peak to below 0.4 keV which is consistent with the data. 
The Fe Ka line is too weak and too broad in this model. 
The same trend can be found when we keep q = 3 fixed 
and leave the outer disk radius as a free parameter. In 
this case, the best-fit (x 2 /d.o.f. = 943/805) results in an 
outer disk radius of only R ont = 4.1^2 r g to model the 
soft excess as smooth as possible. However, we decide to 
fix the outer radius at i? ut = 400 r g and leave q as a free 
parameter. The reflection component represents 21% of 
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Fig. 2.— XMM-Newton EPIC spectra of RBS 320 in the observed 
frame. The 1—12 keV energy range is fitted with a 1PL model. The 
6.4 keV Fe Ka line will be shifted to 4.3 keV at the object's redshift 
of z = 0.495. 



the total 0.2-12 keV output. 

As in the case for RBS 1055, we find a second best-fit 
model (REFL2) with a local fit minimum at an ionization 
parameter of £ = 15001 2 qq erg cm s _1 and q as a free pa- 
rameter, where the reflection component represents 56% 
of the total 0.2-12 keV output. However, this fit mini- 
mum is worse (x 2 /d.o.f. = 981/805) than the fit with the 
low ionization parameter (REFL1) because it contains a 
smaller contribution at hard energies and leaves positive 
residuals around 6-7 keV. 

Applying REFL1G and REFL2G by adding a narrow 
Gaussian component with energy centroid fixed at E — 
6.4 keV and with the line width fixed at a — 0.001 keV 
improves each fit further (Table 5). Leaving the line 
energy as a free parameter only marginally improves the 
fit. The line has EW= 0.05+^3 keV in REFL1G (EW= 

0.06t°;o3 keV for REFL2G). 

The power-law fit of the RGS data (x 2 /d.o.f. = 
220/251) improves when we add two Gaussian line pro- 
files in emission at rest-frame energies oi E — 1.04 ±0.01 
keV and E = 1.24±g;^ keV (x 2 /d.o.f. = 208/246). 
Both lines are unresolved and have rest frame equivalent 
widths of EW= 10t 2 eV and EW= 7tg eV, respectively. 
Because the fit improves by A\ 2 = —8 and -5, respec- 
tively, the lines are detected at the 2-3cr confidence level. 
In the case of the first line, the energy centroid is not far 
from the expected line energy for Ne X, 1.022 keV, but 
thanks to the resolution of the RGS, this energy is ruled 
out at more than 95% confidence. The origins of these 
tentatively detected lines thus remain unclear. Adding 
Gaussian components in emission or absorption at the 
rest-frame energies of expected H or He-like metals re- 
sult in upper limits to the absolute value of equivalent 
widths up to 2 eV. Both disk reflection scenarios fit the 
RGS data equally well. 

5.3. X-ray Spectra of RBS 1124 

A 2PL model again represented the data better 
(x 2 /d.o.f. = 1525/1291) than a 1PL model (x 2 /d.o.f. 
= 1660/1293; see Table 4 and Figure 3). The PLBB 
model does not provide a better fit (k B T = 0.14 ± 0.01 
keV, T = 1.88 ± 0.02, x 2 /d.o.f. = 1580/1291). A partial 
covering model with a column density of Ah ~4x 10 23 
cm~ 2 and a covering fraction of ~21% describes the data 
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very well (x 2 /d.o.f. = 1515/1291). A full-covering and 
partial-covering absorption model with ionized material 
can be excluded because they both require very high col- 
umn densities and ionization parameters (JVh ~ 1 x 10 
cm~ 2 , log (£/(erg cm s _1 ))= 5; see explanation for RBS 
320) . Furthermore, we exclude an Fe edge component at 
7.1 keV (rest frame) to account for the positive residuals 
around 6 keV (rest frame) in 2PL model (90% confidence 
upper limit on r = 0.02). 

Adding a Gaussian component to a 2PL model im- 
proved the fit further by Ax 2 = -9 (2PLG, x 2 /d.o.f. = 
1516/1289). Consequently, the line is detected at the 3a 
confidence level. Fixing the width of the line to a — 0.001 
keV yields x 2 /d-0.f. = 1521/1290. The line is resolved 
since the lower 2a limit of the line width is not consistent 
with the expected value for an unresolved line. Leaving 
the line energy as a free parameter (2PLG2) only im- 
proves the fit marginally (Ax 2 = — 1). 

The self-consistent disk reflection model (REFL2) re- 
sults in a very good fit of the data (emissivity index fixed 
to 3, x 2 /d.o.f. = 1516/1289). We find an extremely high 

ionization parameter of £ = 9000^2000* er g cm s_1 an d 
almost solar abundances. The reflection component ac- 
counts for 75% of the total 0.2-12 keV output. With the 
best-fit parameters of this model, the shape of the reflec- 
tion component is very similar to a simple power-law fit 
except that it models a very broad and smeared out Fe 
line. 

As for RBS 1055 and RBS 320, a second stable local 
X -minimum is found for a fit with low £ values (REFL1, 
X 2 /d.o.f. = 1542/1289). REFL1 introduces discrete line 
features to model the soft excess and the Fe line. Since 
RBS 1124's soft excess is again very smooth, the reflec- 
tion component can only account for 6% of the total 0.2- 
12 keV output in order to avoid introducing significant 
structure to the soft excess. Furthermore, REFL1 pre- 
dicts a slope which is too steep at high energies compared 
to the data. 

Thawing the emissivity index q does not affect the 
quality of REFL1 and REFL2 significantly and q is found 
to be unconstrained in both fit models. Therefore, wc de- 
cide to leave it fixed to q = 3. Considering an additional 
narrow Gaussian line profile (REFL1G and REFL2G) 
with a fixed rest-frame energy of 6.4 keV improves both 
disk reflection fit models marginally. The line has an 
EW= 0.03±g;g| keV in REFL1G and REFL2G. 

After fitting the RGS data with a single power-law 
model (x 2 /d.o.f. = 349/349), no narrow striking spec- 
tral features are visible (upper limits to absolute values 
of EWs on prominent narrow emission and absorption 
lines are up to EW ~ 5 eV). The RGS data are again 
equally well fitted by both reflection scenarios. 

5.4. Joint Fit Analysis of XMM-Newton and Suzaku 

Spectra 

RBS 1124 was also observed with Suzaku on April 14- 
16, 2007 (ObsID 702114010) for a duration of 138 ks (be- 
fore screening for Earth occultation, etc.). The analysis 
and results of this observation are published in Miniutti 
et al. (2010). We include in this paper joint spectral fit- 
ting of the XMM-Newton and Suzaku spectra to achieve 
better constraints on spectral fit parameters for the Fe 
line compared to XMM-Newton alone. In addition, the 
two observations of RBS 1124 gives us the opportunity 




observed energy (keV) 

Fig. 3.— XMM-Newton EPIC spectra of RBS 1124 in observer 
frame. A 1PL is fitted to the 1-12 keV energy band. The 6.4 keV 
Fe Ka line will be shifted to 5.3 keV at the object's redshift of 
z = 0.209. 



to study spectral variability in one of the most luminous 
radio-quiet QSOs on a 13 month time scale and with 
comparable S/N spectra. 

5.4.1. Suzaku RBS 1124 Data Reduction 

The data were obtained with the X-ray imaging CCDs 
(XIS; Koyama et al. 2007) on board Suzaku, which con- 
sist of two front-illuminated (FI) CCD and one back- 
illuminated (BI) CCD (XIS2 was not used), as well as 
the non- imaging hard X-ray detector (HXD; Takahashi 
et al. 2007). RBS 1124 was placed at the "HXD nom- 
inal" position. The XIS and HXD data are processed 
with Suzaku pipeline processing version 2.0.6.13. The 
analysis uses the HEADAS version 6.6.1. software. The 
XIS data reduction follow the guidelines provided in the 
Suzaku Data Reduction Guide 8 . As per a notice for a 
calibration database update 9 the task xispi is run on 
the unfiltered events again. The events are filtered us- 
ing standard extraction criteria (e.g., avoiding the South 
Atlantic Anomaly, low Earth elevation angles, etc.). The 
XIS good integration time is 86 ks per CCD. We extract 
for each of the three XIS CCDs the source events by 
choosing a region 3' in radius, the background events by 
four regions lf5 in radius, and the 55 Fe calibration source 
events. RBS 1124 is detected with a net count rate of 
0.29 ct s- 1 in the XIS FI and 0.39 ct s" 1 in the XIS BI. 
The response (RMF) and auxiliary (ARF) file are pro- 
duced by using the tasks xisrmf gen and xissimarf gen. 
The data for both XIS FI detectors are co-added to form 
a single spectrum using mathpha, addrmf , and addarf . 
Finally the data are grouped to at least 20 counts per 
bin. 

The Mn I Kai and K«2 lines in the calibration spec- 
trum are fitted by two Gaussian components, where the 
K«2 Gaussian component has half the intensity of K«i. 
The nominal rest-frame energies of Mn I Kai and Kct2 
are 5.899 keV and 5.888 keV, respectively, with an ex- 
pected intensity ratio of 2:1. The observed Ka\ line 
peaks at E = 5.897 ± 0.001 keV for the added XIS FI 
and E = 5.897 ±0.002 keV for the XIS BI with an upper 
limit up to a = 16 eV (co-added FI and BI). The FWHM 
energy resolution is 149 eV for the present observation. 

8 http: / /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/abc. html 

9 http: / /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/sci_gain_updatc. html 
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The 12-76 keV HXD/PIN data are also reduced follow- 
ing the Suzaku guidelines. As in Miniutti ct al. (2010), we 
find that an extrapolation of the XIS spectral fit to higher- 
energies does not agree with the PIN data. Walton et al. 
(2010) explain similar Suzaku spectra (e.g., 1H 0419-577) 
with two reflection disk models: one (with kdblur) mod- 
eling the effects of strong gravity in the innermost region 
of the disk, the other (without kdblur) accounting for 
reflection at far distances. Although modeling the data 
better than a 2PL fit, the proposed models, including 
splitting the disk into an inner and outer region, still 
leave strong positive residuals at energies below 20 keV 
(observed frame). Since the origin of the observed steep 
excess in the PIN data below 20 keV (observed frame) 
still remains unclear (likely an additional variable hard 
X-ray source in the HXD/PIN field of view or instrumen- 
tal issues), we decide not to use the PIN data. 

5.4.2. Suzaku X-ray Spectra of RBS 1124 

Miniutti et al. (2010) focus on fitting the Suzaku RBS 
1124 data with the self-consistent disk reflection model 
ref lion and a partially covering model. We briefly sum- 
marize our fit results which also includes detailed power- 
law fits on the Suzaku 0.3-10 keV spectrum of RBS 1124. 
We fit the spectrum over 0.3-9 keV for the XIS BI and 
0.4-11.5 keV for the co-added XIS FI. Furthermore, we 
ignored energies in the range 1.75-1.90 keV to avoid 
calibration uncertainties associated with the instrumen- 
tal Si K-edge. Leaving the normalization difference be- 
tween the XIS FI and BI as a free parameter to account 
for instrumental cross-calibration improves all fits (e.g., 
1PL by Ax 2 = —15). We find a normalization constant 
between both instrument types of approximately 0.98. 

A 2PL model fits also the Suzaku data better than a 
1PL model (Table 4, x 2 /d.o.f. = 958/924 versus x 2 /d.o.f. 
= 1122/926). The PLBB model is a very good represen- 
tation of the data (x 2 /d.o.f. = 947/924, T = 1.72 ±0.02, 
fceT = 0.14 ± 0.01 keV). A partial-covering absorption 
model with neutral material does not explain the data 
well (x 2 /d.o.f. = 968/924, T = 1.94±0.03, A H - 2 x 10 23 
cm~ 2 , and a covering fraction of ^30%). There are no 
obvious narrow discrete absorption features seen that 
would suggest the presence of an ionized absorber as well. 
We return to the 2PL fit and add a Gaussian line profile. 
The 2PLG model improves the fit (x 2 /d.o.f. = 942/921) 
and an Fe line is detected at greater than the 4cr con- 
fidence level. Fixing the line width to a = 0.001 keV 
results in x 2 /d.o.f.= 944/922. The line is not resolved. 

REFL1, the self-consistent disk reflection model with- 
out an additional narrow Fe line, is a very good represen- 
tation of the data (x 2 /d.o.f. = 934/921, emissivity index 
frozen) when an ionization parameter of £disk = 38 
erg cm s _1 is used. The reflection component accounts 
for 13% of the total 0.2-12 keV output. We also find 
in the Suzaku data a second fit minimum incorporat- 
ing a high ionization parameter (REFL2, x 2 /d.o.f. = 
967/922, £ disk = 2200tgoo° er S cm s_1 and keeping the 
emissivity index fixed to 3). REFL2 significantly im- 
proves (A% 2 = -22) when the emissivity index is allowed 
to vary (q = 10t°*, £ disk = 6700t 33 °°* erg cms- 1 ). This 
reflection model makes up 54% of the total 0.2-12 keV 
output. 

We then added a narrow Fe line to the reflection models 
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Fig. 4. — To illustrate the Fe Ka fluorescence line in RBS 1124 
we show the data/model residuals for a 2PL model fit, i.e., without 
the Fe line modeled, to the XMM-Newton EPIC (black - pn, red - 
MOS1, and green - MOS2) and Suzaku (dark blue - XIS FI, light 
blue - XIS BI) spectra. 

(REFL1G and REFL2G) and find improvements to the 
fit compared to REFL1 and REFL2, respectively (Ta- 
ble 5). Model REFL2G improves by Ax 2 = -24 when 
the emissivity index is left as a free parameter. Our fit 
results for model REFL1G are consistent with the val- 
ues given in Miniutti et al. (2010) considering the un- 
certainties. Based on model REFL1G we derive intrinsic 
/o.s-2 = 3.0 x 10- 12 erg cm" 2 s" 1 (/ 2 _ 10 = 4.8 x 10~ 12 
erg cm 2 s 1 ) and a corresponding Lq.5-2 = 3.8 x 10 44 
erg s _1 (I/2-io = 5.8 x 10 44 erg s -1 ). 

The total 0.5-2 keV flux and total 2-10 keV flux both 
increased between the Suzaku and XMM-Newton obser- 
vation, by roughly 60% and 35%, respectively. In the 
context of the 2PL model, there is evidence for both the 
1 keV normalization and photon index of both power- 
law components (soft and hard) to vary between the ob- 
servations. The soft excess contribution according to a 
phenomenological blackbody model (computing the con- 
tribution of the blackbody model to the total 0.5-2 keV 
flux) has decreased from (9 ± 1)% in the Suzaku obser- 
vation to (5 ± 1)% in the XMM-Newton observation 13 
months later. Because the soft excess is observed to vary, 
an origin in a very extended or distant region (greater 
than 13 light months from the central X-ray source) is 
ruled out. The detailed relation between the origins and 
interactions of the soft excess and the hard power law 
for RBS 1124 remain unclear. However, the parameters 
for the Fe line profile do not vary between the two ob- 
servations (line parameters are consistent within the la 
uncertainties), and so we can do a joint fit to achieve 
better constraints on the Fe line. 

5.4.3. Fe Line Constraints Through a Joint Fit 

For the joint fit we use the same energy ranges that 
we used for the single spectral fits. We apply the 2PLG 
model to the data. For the Suzaku data we keep the 
BI/FI normalization constant frozen at c = 0.98. Ty- 
ing the slopes of the soft and hard power laws for both 
observation and only using different normalizations yield 
a very bad fit with very significant residuals. Conse- 
quently, we leave both the 1 keV normalization and the 
photon indices of the soft and hard power law free to vary 
between the Suzaku and XMM-Newton observations. Al- 
lowing the line width and the line energy to vary between 
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Fig. 5.— XMM-Newton EPIC spectra of z = 0.545 RBS 1883 in 
observer frame. A 1PL is fitted to the 0.5-12 keV energy band. 



Table 6. Magnitudes, Rest-frame Luminosity Densities, a, 



Object AB Magnitude/ lo S z 2 500 A log ^ 2 keV a ' 

Name Filter 



RBS 1055 


18.18±0.02/UVVvT* 
17.26±0.03/UVM2 


30.11 


27.41 


1.04 


RBS 320 


15.44±0.01/UVM2 


31.25 


27.15 


1.26 


RBS 1124 


16.44±0.01/UVW1 


30.07 


26.94 


1.20 


RBS 1883 


16.33±0.01/UVW1 
16.29±0.02/UVW2* 










30.98 


26.91 


1.57 



The magnitudes with the superscript "*" are less affected by 
strong UV emission lines and are used for the computation of the 
rest- frame luminosity at 2500 A and a x- The effective wavelengths 
for the filters are: UVW2 - 2120 A, UVM2 - 2310 A, and UVW1 - 
2910 A. The rest-frame luminosity densities ' 2500 £ and I2 keV arc 

given in units of erg s Hz - 1 . 



both observations does not improve the fit further. 

We derive an Fe line width of (T FGjl . ost = 0.141°; Jg keV 
from the joint fit (Figure 4). If we omit the line com- 
ponent, the fit worsens by A\ 2 = 24, while fixing the 
line width to 1 eV worsens the fit by A\ 2 = 8. Conse- 
quently, using a joint fit to both observations, we detect 
the Fe line with a much higher significance than in the 
individual observations. The line is resolved and has an 
intensity of I = (5.2 ± 2.1) x 1(T 6 ph cm" 2 s _1 . 

5.5. X-ray Spectra of RBS 1883 

The observation of RBS 1883 suffers heavily from very 
high background count rates throughout almost the en- 
tire observation. Although the low number of counts does 
not allow us to study RBS 1883 as detailed as the other 
RBS QSOs presented in this paper, a 2PL model fits the 
data much better than a 1PL (x 2 /d.o.f. = 180/86 versus 
119/84; see Table 4 and Figure 5). A PLBB model re- 
sults in a fair fit of the data {k B T = 0.13±0.01; x 2 /d-cf. 
= 124/84). 

A partial covering model improves the fit (x 2 /d.o.f. 
= 113/84) by finding a very high covering fraction of 
(82±6)% and A H , CO v = (l-C^) x 10 23 cm" 2 . The S/N 
of the data does not allow us to test for a warm absorber 
in RBS 1883; the ionization parameter and iVH, W arm stay 
unconstrained. For the same reason we are not able to 
apply a reflecting disk model to the data. 

Although no emission line features, e.g., an Fe Ka line, 
are visible in the spectra, we included a Gaussian compo- 
nent fixed to an energy of 6.4 keV (rest frame) to give up- 
per limits on the Fe Ka emission line. For an unresolved 
Fe line (width a = 0.001 keV), we find an 90% confidence 
level upper limit on the EW< 0.66 keV (EW< 0.84 keV 
for a resolved Fe line with a fixed to 0.4 keV). 

6. OPTICAL-TO-X-RAY SPECTRAL INDICES 

The optical monitor (OM) onboard XMM-Newton pro- 
vides UV photometric data for all objects taken simul- 
taneously with the X-ray data and therefore enables 
us to calculate optical-to- X-ray spectral energy distri- 
butions (a ox ) without the uncertainties introduced by 
source variability. The OM source lists produced by 
the XMM-SAS pipeline contain calibrated source mag- 
nitudes in the AB system that are also corrected for the 
long term degradation of the OM. In addition, we de- 
redden the UV data using the maps by Schlegel el al. 



(1998) and the wavelength dependencies by Cardelli et 
al. (1989). 

The four OM observations of the RBS-QSOs were 
taken with different filter combinations. To derive the 
fc-corrected rest-frame luminosity densities at 2500 A, 
we use the mean UV continuum slope a uv = —0.44 (fre- 
quency index) of the SDSS QSO spectra (Vanden Berk 
et al. 2001). If we have more than one magnitude mea- 
surement for an individual object, we select the filter 
observation which is less affected by a strong UV emis- 
sion line. The only OM detection of RBS 320 in filter 
UVM2 coincides with the C IV line in the object and the 
derived values should be considered with care. The rest- 
frame luminosity densities at 2 keV is computed from the 
0.5-2 keV rest-frame luminosities given in Table 3 and 
the photon indices measured in the same energy range 
(Table 8). Considering the redshift of our objects, the 2 
keV rest-frame luminosity density will be redshifted into 
the 0.5-2 keV band. The measured magnitudes and the 
derived quantities are listed in Table 6. 

The derived a ox values cover the same range that is 
found by Grupe et al. (2010) in a sample of soft X-ray se- 
lected AGN drawn from the RASS survey and observed 
with Swift. They also observed RBS 1124 and find an 
a ox = 1.23. Considering this Swift sample, RBS 1055 
yields an a ox value on the low edge of the distribution. 
Our values are also consistent with the sample by An- 
derson et al. (2007, 2003) who spectroscopically identi- 
fied 7000 RASS AGN from the Sloan Digital Sky Survey 
(data release 5). 

7. LONG-TERM FLUX VARIABILITY OF THE RBS-QSOS 

Since all the studied objects in this work have been 
observed with ROSAT and XMM-Newton separated by 
almost two decades, both measurements allow a crude 
long-term flux variability study. The RBS contains the 
flux-brightest (with some spatial restrictions) objects 
from the RASS which was observed in 1990/1991 dur- 
ing the first half year of the ROSAT mission. RBS 1124 
has also been observed by the XMM-Newton Slew Sur- 
vey, Suzaku, and Swift (see Miniutti et al. 2010). RBS 
320 was observed four years after the XMM-Newton ob- 
servation with Swift on 2008 November 7 (target name: 
QSO B0226-4110; target ID: 37516). We used the online 



12 



Krumpe et al. 



6 












; A - - . 


- - _ ^RBS 


1124 




' m 5 








~ ~ - □ ; 










l ' ! 


OB 

Sh 
CO 

2 3 

o 


: a — 




1055 


:/ ; 
& i 

------- -i 


f 

0.5-2 keV 


i a - 


RBS 

RBS 


320 

1883 


□ o ; 














1990 1995 2000 2005 2010 

Year 

Fig. 6. — Intrinsic (Galactic-absorption corrected) flux vs. 
time for the studied RBS-QSOs. The symbols represent measure- 
ments from different missions: triangle - ROSAT, square - XMM- 
Newton, diamond - Swift, asterisk - XMM-Newton Slew Survey, 
and cross - Suzaku. 



Swift XRT data reduction pipeline 10 (Evans et al. 2009) 
to extract a ~13 ks spectrum. Binning the data to 20 
counts per bin, the low S/N (^1000 counts) Swift spec- 
trum requires only a single power-law fit (r = 2.58±0.09, 
X 2 /d.o.f. = 41/42) with no intrinsic absorption. This is 
consistent with the 1PL fit of the XMM-Newton data 
considering the uncertainties. 

The 0.5-2 keV ROSAT fluxes are taken from the RBS 
catalog (Schwope et al. 2000). They assumed a photon 
index of T = 2 to compute count-to-energy conversion 
factor for the ROSAT PSPC detectors considering the 
Galactic column density toward the X-ray source. To 
keep possible flux changes due to the different flux com- 
putation methods between R OSA T (RBS) and the mea- 
surements from other missions as low as possible (flux 
differences up to 10%), we applied the RBS technique to 
all observations. Figure 6 shows all measured fluxes for 
the studied RBS-QSOs versus time. Except for the de- 
crease and rise of RBS's 1124 flux between 2006 Novem- 
ber and 2008 May, all objects show only flux variability 
on a level of ~5%-20% based on the few data points 
available. 

8. CONSTRUCTION OF AN AVERAGE X-RAY SPECTRUM 

We decided to construct an average X-ray spectrum 
for the four sources studied here. First, the average spec- 
trum may unveil a possible relativistic component of the 
Fe Ka emission line in the most luminous radio-quiet 
RBS-QSOs due to the increase in S/N. Secondly, the av- 
erage broadband spectral properties in these objects can 
be compared to the average broadband spectral prop- 
erties obtained for high-redshift, low S/N spectra with 
very similar X-ray luminosities found routinely by the 
XMM-Newton and Chandra surveys (e.g., Mateos et al. 
2010; Brusa et al. 2005). This will test if our most lumi- 
nous radio-quiet RBS-QSO X-ray spectra, which allow 
us to determine their spectral properties directly, are the 
standard cases of luminous radio-quiet QSOs at high red- 
shifts. 

1(1 http : //www . swift . ac . uk/user_obj ect s/ 



The averaging process follows Corral et al. (2008) with 
some differences. The original process was designed to 
average a large number of sources to improve the S/N, 
whereas here we are only averaging four sources. There- 
fore, the influences of individual sources on the average 
spectrum is considerably higher. We tested the source- 
to-source influence on the average spectrum by generat- 
ing sets of average spectra, removing one source in each 
case. The results are robust, i.e., the derived properties 
vary within their 2a limits. In particular, removing RBS 
1055, which has the broadest Fe line, does not change 
the properties of the averaged Fe line. 

We are forced to exclude the pn data for RBS 1883 
since, given the small number of counts at high energies 
in this case, it distorts the average spectrum shape at 
those energies. The main steps of the averaging process 
are as follows: We fit a 2PL model to each spectrum, 
instead of the simple absorbed power law as described in 
Corral et al. (2008), since, as we show above, this is a 
very good representation of the X-ray data for the whole 
sample. Then we obtain the incident spectrum for each 
source, i.e., the source spectrum before entering the de- 
tector, by using the previously obtained model and the 
Xspec command euf spec. Doing so, we derive the source 
spectra in flux units of keV cm~ 2 s _1 keV -1 . Next, we 
de-absorb each spectrum for the Galactic absorption at 
each source position and shift them to their rest frame. 
To combine the individual spectra, we rescale them by 
forcing the rest frame 2-5 keV flux to be the same for 
all objects. Finally, we rebin the spectra to a common 
energy grid in a range of 1-13 keV, ensuring that all 
bins contain at least 600 real counts, and average them 
by using a standard mean. Individual errors on the real 
spectra are propagated as Gaussian along the whole av- 
eraging process. The significances of the individual spec- 
tral features in the final average spectrum are quantified 
by using simulations. We simulate each source 100 times 
using its best-fit model and keeping the same spectral 
quality and flux as for the real source. We apply exactly 
the same averaging process to the simulated spectra as we 
do for the real spectra. From the 100 simulated spectra, 
we compute the confidence limits (la and 2a confidence 
encompass 68% and 95% of the simulated values, respec- 
tively) and a simulated continuum, which account for the 
mixture of 2PL shapes. As already suggested by the in- 
dividual X-ray spectra, the resulting average spectrum 
contains in comparison to a 2PL no spectral features ex- 
cept an Fe line which is detected above the 2a confidence 
limit. The Fe K bandpass is shown in Figure 7. 

Following the analysis on individual objects, we first 
fit a 2PL to the 1-12 keV rest frame data. We achieve 
a good fit of x 2 /d-o.f - 22/19 (r soft - 3.4±g'J, r hard = 

1.7^ ' 2 )- The remaining positive residuals can be well 
fitted by adding a Gaussian profile. The 2PLG model 
improves the fit to x 2 /d.o.f. = 8/16 and yields r so f t = 
3.4±g:i> A soft = 0.8±%% ph cm^ 2 s" 1 keV" 1 (normaliza- 
tion at 1 keV), r hard = 1.7lg;J, A hald = 0.7t°; 2 ph cm" 2 
s^ 1 keV~\ £ F c,rcst = 6.4 ± 0.1 keV , o- F c,rcst = 0.2 ± 0.2 
keV, and J Fe = 6.004±0.002 ph cm~ 2 s" 1 .' We find an av- 
eraged rest-frame equivalent width of EW= 0.12 ± 0.06 
keV in the stacked spectrum of the RBS-QSOs. The 
fit still leaves some tempting looking positive residuals 
around 4.5-6 keV (see Figure 7), but adding an addi- 
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narrow fluorescent Fe Ka line. Page et al. (2004b) study 
the spectral properties and the spectral energy distri- 
bution of five high-luminosity radio-quiet QSOs. Their 
sample spans a redshift range of z = 0.24 — 2.04. The 
two objects with a similar redshift range to our objects 
have lower 2-10 keV luminosities than our RBS-QSOs. 
Piconcelli et al. (2005) investigate the X-ray properties 
of a much larger sample of 40 PG quasars in a red- 
shift range of 0.036-1.718. The sample contains only 
five radio-loud QSOs and covers an X-ray luminosity of 
10 43 < io.s-io/Icrgs- 1 ] < 5 x 10 45 . Therefore, PG 
quasars share very similar properties to our sample and 
we will compare our results mainly to this sample. 

None of our studied RBS-QSOs require the presence 
of a full covering cold absorber along the line of sight in 
excess of the Galactic equivalent hydrogen column den- 
sity. This is consistent with the findings of Piconcelli et 
al. (2005) where only 3/35 radio-quiet PG quasars show 
evidence for the presence of a cold absorber. The optical 
spectra of our RBS-QSOs reveal AGN-typical features: 
broad emission lines and UV-excess. These type of ob- 
jects (Seyfert I classification) are, in general, known to 
agree well in optical and X-ray properties with the classi- 
cal AGN unification scheme (Antonucci 1993). However, 
RO SAT's sensitivity to the soft energy band (<2.4 keV) 
Table 7. 1PL Photon Indices for. Different Redshifts of the Averaged selects against X-ray absorbed AGNs. Therefore, our 



1.30 



1.20 



1.10 



1.00 



0.90 



0.80 




Energy (keV) 

Fig. 7.— Added XMM-Newton spectrum of RBS 1055, RBS 320, 
RBS 1124, and RBS 1883. The stack contains ~35,000 counts in 
the shown 3-9 keV range. The squares represent the ratio between 
the flux of the stacked spectrum and the simulated continuum spec- 
trum. The uncertainties of the bins of the stacked spectrum are 
shown as error bars. The dotted lines show the la confidence lim- 
its, while the dashed line represents the 2cr confidence limits. The 
solid gray line at a ratio of one illustrates the continuum flux. 



RBS-QSO Spectrum 



z 


Photon index 


X 2 /d-o.f. 




(1-12 keV) 


0.0 


2.05±0.02 


167/21 


0.5 


1.91±0.03 


45/22 


1.0 


1.85±0.04 


22/21 


1.5 


1.82±0.04 


20/20 


2.0 


1.81±0.05 


20/19 



tional Gaussian profile docs not yield a statistically sig- 
nificant improvement in the fit. In order to give at least 
an upper limit for the equivalent width of a broadened 
disk line, we add a laor profile to the 2PLG model. We 
only allow the normalization of the laor disk line to vary 
(i? in fixed to 1.235 r g ) and find EW< 340 eV. 

Many other studies stack X-ray spectra of high- redshift 
AGNs and yield a data quality which is sufficient only to 
derive the average broadband spectral properties. To 
compare our data to these studies, we fit a 1PL model 
only to our stacked RBS-QSO spectrum. Moreover, we 
derive the observed power-law indices in the 1-12 keV 
range depending on different redshifts in Table 7. For 
z < 0.5, the 1PL fit is a bad representation of the data, 
but once the soft excess is shifted outside the 1-12 keV 
range this model results in a good fit. Above z ~ 2 the 
Compton-hump is expected to be redshiftcd into the 1-12 
keV range. Since the Compton hump strengths for each 
of our objects remain unclear, we can only speculate that 
the added spectrum may show further hardening. 

9. DISCUSSION 

Although we investigate only a small sample, all of 
our objects show very homogenous spectral properties: 
a very smooth soft excess and a hard power-law com- 
ponent. Furthermore, in the three high S/N spectra of 
RBS 1055, RBS 320, and RBS 1124 we detect a weak 



sample is biased toward unabsorbed AGNs. 

Ionized absorbers are found routinely in ^50% of 
Seyfert galaxies (Reynolds 1997). Piconcelli et al. (2005) 
detect absorption features due to ionized gas in 18 out 
of their 40 PG quasars with similar and lower X-ray lu- 
minosities. They conclude that the detection rate of ion- 
ized absorbers appears to be independent of X-ray lu- 
minosity. For each of our objects we find no evidence 
for the presence of ionized absorbers. The RGS data of 
the three high S/N RBS-QSOs also show no evidence for 
ionized absorbers. The S/N of our spectra is comparable 
to Piconcelli et al. (2005). Other examples of ionized ab- 
sorbers in luminous radio-quiet AGNs are found by Page 
et al. (2004b) (one out of five objects at z = 0.24) and 
Reeves et al. (2003) (z = 0.18). The extreme smoothness 
of our RBS-QSO spectra in the soft energy band, with no 
signs for ionized absorbers, is not consistent with these 
previous studies and may suggest a deficiency of ionized 
absorbers in the very high end of the QSO luminosity 
function. 

Partially covering models are only able to sufficiently 
describe the observed data of RBS 1124 and RBS 1883. 



However, in both cases this requires 7Vh. ( 



10 



2.3 



which results in very high bolometric luminosity (a few 
10 46 erg s -1 assuming the bolometric correction to the 2- 
10 keV intrinsic luminosities of Marconi et al. 2004) and 
consequently in supcr-Eddington ratios of Lboi/-^edd=5.3 
for RBS 1124 and 3.9 for RBS 1883. Therefore, a par- 
tially covering model seem to be unlikely to explain the 
X-ray spectra for these two objects. 

9.1. Soft Excess 

The need for a 2PL model to fit the overall X-ray spec- 
tra well in all objects indicates that all objects have a 
significant soft excess. Moreover, in all objects in the 
current sample, the soft excess is well fit by a soft power 
law in CCD detectors. The RGS spectra of the three 
highest S/N spectra are also fit well by a simple power 
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law, with no strong evidence for narrow emission or ab- 
sorption lines. The vast majority of luminous AGNs is 
known to show a soft excess (e.g., Page et al. 2004b). 
For comparison, we derive the photon indices by select- 
ing only the energy ranges of 0.5-2 keV and the 2-12 keV, 
respectively. We list the results of a single power-law fit 
to these energy ranges in Table 8, i.e., these indices are 
not based on our 1PL or 2PL model fits. Our mean soft 
photon index ((r so ft) = 2.36) is harder than the aver- 
age PG quasar value of (r soft ) = 2.73. Only RBS 320 
and RBS 1883 have soft photon indices comparable to 
the average PG quasar value, but these objects also ex- 
hibit much softer 2-12 keV photon indices than the other 
sources. 

Miniutti et al. (2009) evaluate the strength of the soft 
excess for the PG quasars based on a phenomenological 
blackbody model by computing the contribution of the 
blackbody component to the 0.5-2 keV total flux finding 
an average of ^30%. We list this quantification of the 
soft excess for our RBS-QSOs in Table 8. Except RBS 
1883, all our objects show a much smaller soft excess 
strength (2%-9%) than the, on average, lower luminous 
PG quasars. 

In the PLBB fits of our objects, we derive a mean 
(& b 7bb) = 0.12 keV with a dispersion of a = 0.02 keV. 
This is consistent with the values found for AGNs with 
intermediate- mass black holes (and -Lo.5-10 ~ 10 43 erg 
s _1 ) from Miniutti et al. (2009) and slightly lower (but 
still consistent) than the mean value for PG quasars 
((/c b Tbb) = 0.136 keV, a = 0.021 keV). As pointed out 
by several authors (e.g., Gicrlihski & Done 2004; Ross 
& Fabian 2005; Crummy et al. 2006), the observation of 
the relatively sharp blackbody temperature distribution 
of ^0.10-0.15 keV over such wide range of X-ray lumi- 
nosities does not support the idea that the soft excess in 
AGNs is connected to the high-energy tail of the thermal 
emission from the accretion disk. We find similar black- 
body temperatures in the PLBB fits to our RBS-QSOs. 
In this context, the soft excess consequently has to be ex- 
plained by a model that produces very similar Tbb values 
independently of AGN X-ray luminosity when modeled 
by a blackbody component, as well as reduces the soft 
excess strength from ^30% for low and medium X-ray 
luminous to ~5% for high X-ray luminous radio-quiet 
QSOs. 

9.2. Hard X-ray Continuum 

The mean 2-12 keV photon index of our RBS-QSOs 
is (Fhard) = 1-96 (<r = 0.30). Taking into account that 
this value is only based on four objects, it agrees well 
with the mean hard power-law component of radio-quiet 
PG quasars ((r hard ) = 1.89, a ~ 0.36) and the QSO 
sample of Page et al. 2004b ((r hard ) = 1.93 ± 0.02). Pi- 
concelli et al. (2005) find some significant (>99% con- 
fidence level) correlations among the 2-12 keV power 
law/luminosity and QSO physical properties. Our ob- 
jects, testing these correlations at the high end of the 
QSO luminosity function, fit well into the Mbh~£2-io 
and FWHM(H/3)-T correlation. They are somewhat 
consistent with the T-Mbh correlation, but do not 
fit well into the FWHM(H/3)-L 2 -io correlation (higher 
FWHM(H/3) values are expected for our luminosities). 
The hard photon indices of our individual objects are 
consistent with the r-Lboi/^edd correlation (e.g., Lu & 



Table 8. Additional Fit Results from the RBS-QSO X-ray Spectra 



Object 


r 


r 


Flux Ratio 




Name 


(0.5-2 keV) 


(2-12 keV) 


(bb/0.5-2 keV) 


(2-10 keV) 


RBS 1055 


1.81±0.02 


1.62±0.03 


n nis+ - 002 
u.oi»_ 004 


24 


RBS 320 


2.69±0.03 


2.13±0.07 


0.05±0.01 


8.7 


RBS 1124 x 


2.00±0.02 


1.73±0.03 


0.05±0.01 


8.0 


RBS 1124 s 


1.93±0.03 


1.69±0.03 


0.09±0.01 


5.8 


RBS 1883 


2.93±0.20 


2 36+ ' 48 

Z.OD_ () 44 


26+ ' 01 


5.5 


Mean 










(XMM only) 


2.36±0.47 


1.96±0.30 




11.6 



Note. — The flux ratio bb/0.5-2 keV quantifies the 
Galactic absorption-corrected contribution of the soft 
excess to the total observed 0.5-2 keV flux; this was 
computed using a PLBB fit. The Galactic absorption- 
corrected 2-10 keV rest-frame luminosity is based on the 
best-fit model of each object (in units of 10 44 erg s _1 ). 
Superscript "X" refers to the XMM-Newton observation 
of RBS 1124, while "S" indicates the Suzaku observa- 
tion. For the error of the mean photon indices we list 
the dispersion a. 



Yu 1999; Shemmer et al. 2008). 

9.3. Origin of the Fe Line 

Three out of our four studied objects have spectra with 
a sufficiently high S/N to detect an Fe line whose ener- 
gies are consistent with a neutral fluorescent Fe Ka line 
at 6.4 keV. A line energy of 6.7 keV from Hc-like Fe 
emission can be excluded in the individual and stacked 
XMM-Newton observations with a statistical significance 
of 1.5(7 to >3<t. All of the lines are either barely resolved 
or nearly resolved (<7F C , r cst = 0.15 — 0.5 keV). The rest- 
frame equivalent widths are 70, 130, and 140 eV (XMM- 
Newton observations only) and are consistent between 
the individual objects taking into account the uncertain- 
ties. The average equivalent width of the individual ob- 
servations of (EW) = 110 eV is in very good agreement 
with the derived equivalent width of the stacked RBS- 
QSO spectrum (EW= 120 ± 60 eV). 

Jimcncz-Bailon et al. (2005) investigate the properties 
of the fluorescent Fe Ka line in the PG quasar sample 
of Piconcelli et al. (2005). A narrow Fe line is found in 
50% of all radio-quiet PG quasars. The vast majority of 
their narrow lines have an energy centroid consistent with 
E = 6.4 keV, indicating matter in low ionization states 
(Fe I-Fe XVII) which is expected to be located at large 
distances from the central X-ray source. However, at lu- 
minosities L2-10 greater than 3 x 10 44 erg s _1 , they only 
have one radio-quiet object with sufficient S/N to detect 
a line. Our study adds three data points in this poorly 
studied luminosity range. The detection of only narrow, 
yet barely resolved, Fe lines from neutral/low ionization 
stages agree well with the results of PG quasars and other 
QSO samples (e.g., Page et al. 2004b; Porquet et al. 
2004). Therefore, narrow Fe emission from low ioniza- 
tion states is a common spectral feature for low to very 
luminous AGNs and suggests that a common mechanism 
or set of mechanisms governing transmission and/or re- 
flection processes must exist across a wide range of lumi- 
nosities. 

We compute the FWHM velocity line width for the Fe 
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Ka lines in each of our objects based on the 2PLG fits 
(Tabic 4). For RBS 1124, we used the line width from the 
joint Suzaku and XMM-Newton fit since it gives slightly 
better constraints and the values from both observations 
agree within their uncertainties. Assuming that the ma- 
terial in which the Fe line originates is in a Keplerian 
orbit, the FWHM velocities of 53, 000 ± 20, 000 km s" 1 
(RBS 1055) and 16,000i^;°o° km s_1 ( RBS 1124 ) cor- 
respond to distances of only 2.4jlf f and 4:.l±l 4 7 9 light 
days, respectively, from the central supermassive black 
hole (using (v 2 ) = |w| WHM , Netzer 1990). These dis- 
tances are even closer to the black hole than the Broad 
Line Region (BLR) that is believed to emit the broad 
optical emission lines. If the velocities are Keplerian the 
Mg II emission in RBS 1055 and R/3 in RBS 1124 origi- 
nate from radii of ^220 and ~58 light days, respectively. 
For RBS 320 the lower 90% confidence limit is r > 1.2 
light days (<75,000 km s -1 , with H/3 emission originat- 
ing at a distance of ^1300 light days). 

It is somewhat surprising to detect neutral Fe emission 
from a radius of only a few light days, hence closer than 
the BLR, in the most luminous QSOs. It is reasonable to 
expect that the material may be highly ionized at these 
radii. Shielding processes could sustain neutral material 
so close to the central engine. However, a partial covering 
model does not fit the data of RBS 1055 well and is also 
unlikely for RBS 1124 (see discussion in Section 9). An- 
other possible scenario of maintaining neutral material 
very close to very luminous AGNs may be non-isotropic 
emission of the ionizing continuum radiation. 

The PG quasar sample shows that the detection of rela- 
tivistically broadened Fe lines is not common for high lu- 
minous AGNs except in three cases. Our non-detections 
are consistent with the notion that relativistically broad- 
ened Fe lines are rare at the high end of the QSO lu- 
minosity function. This finding seems to also apply to 
high luminosity type II QSOs for which Krumpe et al. 
(2008) do not detect a relativistically broadened Fe line 
(EW upper limit of 1 kcV) in the stacked spectrum of 
type II QSOs, while lower luminosity AGNs do show ev- 
idence for such a line (EW= 2.0+g 7)- Possible selection 
effects due to the limited S /N in the spectra of the lumi- 
nous QSOs may hamper the detection of broad Fe lines. 
We determine the 90% confidence level upper limit of a 
laor component for each of RBS 1055, RBS 320, and 
RBS 1124 and find that broad Fe lines with EWs up to 
~ 140 — 400 kcV could be hidden in the spectra. The 
three detected broad Fe lines in the PG quasar sample 
have EW= 100-500 eV (in objects with L 2 -w < 3x 10 44 
erg s" 1 ). De la Calle Perez et al. (2010) investigate a 
sample of radio-quiet AGNs and QSOs and detect only 
for 2 QSO (<6% of all QSOs) a relativistically broad- 
ened Fe line. They also conclude that the measured 
EW is always below 300 eV and has a mean value of 
(EW) = 143 ± 27 eV. An extreme case is the radio-quiet 
very luminous QSO PG 1247+267 (Page et al. 2004b, 
z = 2.04, L 2 -w = 9 x 10 45 erg s _1 ) which shows an 
apparently broad line with EW= 421 ± 215 eV. There- 
fore, the expected EWs of broad Fe lines based on pre- 
viously published results are not detectable in our RBS- 
QSO spectra. At low redshifts, low luminosity Seyfert 
AGNs show a relativistically broadened Fe line in ~50% 
of the objects (Nandra et al. 2007). However, the major 



fraction of the detected lines in these high S/N spectra 
has EW< 200 eV and therefore falls below our detec- 
tion limit for a broad Fe line. Consequently, our results 
do not allow a conclusion on similar or different physical 
properties in the accretion disk very close to the central 
engine in low and high luminous AGNs. 

9.4. Comparison to Low-z AGNs 

Low- z AGNs vary significantly in their individual X- 
ray spectral components and feature a wide range of fre- 
quency of partial coverers, warm absorbers, soft excess 
shapes and strengths, and/or Compton reflection hump 
strengths. It is therefore not straightforward to compare 
our results to a somewhat typical low-z AGN to evaluate 
what changes have to be applied to the X-ray spectrum 
(except the much higher luminosity) to turn a low-z AGN 
into one of the most luminous QSOs. 

Since we are left with the situation where our RBS- 
QSO spectra can be well fitted with power law and self- 
consistent reflection disk models (EPIC and RGS data), 
we decided to compare our findings with Crummy et al. 
(2006). They used a sample of 34 mainly low-redshift 
Seyfert I galaxies and type I AGNs to show that in 25 
out of 34 sources, relativistically blurred disk reflection 
models are a significantly better fit than a PLBB model. 
All three of our high S/N RBS-QSO observations are bet- 
ter fitted with reflection disk models than PLBB mod- 
els. However, a 2PL fit is almost as good as a reflection 
model, because the extremely smooth soft excess in the 
EPIC and RGS data is well fitted by a soft power law. A 
reflection model has to accomplish the modeling of the 
smooth soft excess similar to a power law and requires a 
high degree of relativistically blurring to smear out the 
sharp soft X-ray emission lines caused by reflection from 
the disk. This is achieved by either very high ioniza- 
tion parameters, small inner disk radii, or high emissiv- 
ity indices (causing the emission of the major fraction of 
the radiation within very small inner radii), or a com- 
bination of those. In this case, the limited parameter 
space no longer allows modeling of the Fe line in a sat- 
isfying way as in REFL1 and REFL2, requiring us to 
include the extra, narrow Fe component as in REFL1G 
and REFL2G. Moreover, we find a degeneracy in the disk 
reflection model between accretion disk states at seem- 
ingly opposite physical extremes: REFL1 has a very low 
ionization parameter and accounts only for ~5%-20% of 
the total 0.2-12 keV flux, while REFL2 exhibits a very 
high disk ionization parameter (~35%-75% of the total 
output). 

Some objects in the Crummy et al. (2006) sample (e.g., 
TON S180, where the soft excess is also better modeled 
with a power law rather than a blackbody model) also 
have a very smooth soft excess and require high ioniza- 
tion parameters and high emissivity indices. Crummy et 
al. (2006) allow the emissivity index to vary freely for 
all their sources and find an index higher than 6 for 22 
sources. Only two of their sources are consistent with an 
index below or equal to 3. Leaving the emissivity index 
free to vary in our fits results in unconstrained values in 
most of the cases due to the limited S/N. Only for RBS 
320 and REFL2 of the RBS 1124 Suzaku observation we 
achieve well constrained emissivity indices higher than 6. 

Crummy et al. (2006) find very small inner disk radii 
indicating fast rotating Kerr black holes which is in 
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agreement with our data. Interestingly, their 0.3-12 keV 
flux fraction for the disk reflection models is always above 
0.25 which contradicts our models REFL1 and REFL1G 
(low ionization parameter; flux fraction ~5%-20%). The 
best-fit ionization parameters in Crummy et al. (2006) 
cover a wide range including the values found for our 
RBS-QSOs. They also find that most sources have sim- 
ilar spectral indices with both models (PLBB and disk 
reflection), while several sources have higher indices with 
the disk reflection model. The ranges of our PLBB and 
REFL1G photon indices agree well with the correlation 
given by Crummy et al. (2006). For REFL2G we find 
lower spectral indices than with a PLBB. 

None of the degenerate disk reflection scenarios, i.e., 
low/high ionization parameter, for our RBS-QSOs can 
be favored or ruled out. The limited S/N of the RGS 
data also does not allow any constraints. If the soft ex- 
cess in these objects is not caused by reflection in the 
innermost disk, there is in general no need for relativis- 
tical disk reflection models in our high luminosity QSOs 
as the observed Fe line does not originate from the in- 
ner disk although we cannot rule out the possibility that 
there may exist a weak broad Fe line. Either the ba- 
sic assumptions of reflection models (e.g., geometrically 
flat accretion disk, gas of constant density, etc.) are not 
valid for the most luminous QSOs, or different sets of 
models are required to describe their physics. If, on the 
other hand, reflection models turn out to be the right 
interpretation for very luminous QSOs, this would sup- 
port the assumption by Page et al. 2004b that high lumi- 
nosity radio-quiet QSQs have identical X-ray continuum 
properties as those found in lower-luminosity AGNs and 
Seyfert I galaxies and therefore similar physical processes 
are likely to govern accretion disks over a wide range of 
luminosities and accretion rates. 

9.5. Comparison to High-z QSOs 

In this section we want to explore if our stacked high 
S/N RBS-QSO spectrum can serve as a template spec- 
trum for low S/N spectra of luminous QSOs at higher 
redshifts (z > 1). The PG quasar sample of Piconcelli et 
al. (2005) contains five objects with z > 1 (up to z = 1.7) 
and X-ray luminosities of £2-10 ~ (15 — 130) x 10 44 erg 
s _1 . These objects have an average 2-12 keV photon in- 
dex (r 2 _i2) = 2.1 with a typical error of ±0.3. Our pho- 
ton index for the stacked RBS-QSO spectrum shifted to 
z = 1.0 is Ti-12 = 1.85 ± 0.04. Taking into account the 
slightly different energy range for the determination of 
the photon indices and the involved uncertainties, both 
values are in good agreement. 

Although we are aware of the large uncertainties in us- 
ing a sample of only three objects, we derive the average 
soft and hard photon indices for the three z > 1 ob- 
jects of Piconcelli et al. (2005). Their (r soft ) = 3.4 and 
(Fhard) = 1-4 agree well with the 2PL fit to our stacked 
RBS-QSO spectrum (r soft = 3.4±g;|, r hard = 1.7jft£). 

The limited S/N for luminous QSOs at high redshifts 
restricts the knowledge that can be gained from these ob- 
jects. By stacking the individual X-ray spectra the S/N 
is improved significantly and the added spectrum can be 
used for further analysis which is not feasible on indi- 
vidual spectra. However, the interpretation of stacked 
spectra a priori assumes that the intrinsic X-ray prop- 



erties in the sample are very similar. Radio-quiet AGNs 
correspond to the vast majority (<~90%) of all AGNs. 
Although limited to only four objects at the present, 
our RBS-QSOs are found to have surprisingly homoge- 
nous properties in their X-ray spectra. This confirms the 
stacking technique assumption of very similar properties 
in the individual spectra at least for very luminous QSOs. 

In one of the largest samples, Mateos et al. (2010) 
stack the spectra of 487 type I AGNs in the XMM- 
Newton Wide Angle Survey (XWAS) in different red- 
shift bins. The two highest redshift bin subsamples with 
(z) = 2.01, 2.56 and L 2 - W - 4 x 10 44 erg s" 1 both yield 
a best-fit photon index of (r) = 1.81 ±0.05. Our stacked 
spectrum redshifted to z = 2 is best described by a single 
power-law fit with T — 1.81±0.05. Although both results 
are in perfect agreement, Mateos et al. (2010) mention 
that their value could be biased in such a way that at 
high redshifts the efficiency of detecting hard objects is 
increased, while for soft sources it decreases. They also 
point out that above ~ 10 44 erg s _1 , the average photon 
index seems to remain fairly constant over two orders of 
magnitude in luminosity. Mateos et al. (2010) find little 
or no evidence for X-ray absorption (only ~3% of the 
sources show significant full covering cold absorber along 
the line of sight) and a soft excess model (PLBB model 
versus single power law) represents the data better in 
^36% of the objects below z = 0.5 (F-test significance 
>99%, note that their median number of total EPIC 
counts in the 0.2-12 keV band in the individual spec- 
tra is ^300 counts). Our objects fit well the correlation 
of the soft excess luminosity versus intrinsic power-law 
continuum luminosity presented by Mateos et al. (2010). 

Corral et al. (2008) stack the X-ray spectra of 606 type 
I AGNs in the XWAS and study the properties of the 
average Fe line. For the highest luminosity sub-sample 
£0.5-2 ~ (6 — 660) x 10 44 erg s _1 , they find an average 
Fe line EW of 50^° eV. Our average Fe line EW from 
the stacked spectrum is EW= 120 ± 60 eV, which is con- 
sistent with their value within the uncertainties. Chaud- 
hary et al. (2010) use 507 AGNs from the 2XMM catalog 
and find for subsamples with L2-10 ^ 5 x 10 44 erg s _1 
in a redshift range of z = 0.8 — 5.0 upper limits on the 
EW also consistent with our values. These new stacking 
results and our direct measurements for the RBS-QSOs 
significantly disagree with the large EWs found by Stre- 
blyanska et al. (2005) and Brusa et al. (2005) using the 
stacking technique. 

9.6. X-ray Baldwin Effect 

A unique and important contribution of our study is 
that we can test the X-ray Baldwin effect, the anti- 
correlation between X-ray luminosity and EW of the nar- 
row Fe line (Iwasawa & Taniguchi 1993), for the most 
luminous QSOs. Several previous studies have lively de- 
bated the existence and the significance of the X-ray 
Baldwin effect. For example, Nandra et al. (1997), Page 
et al. (2004a), and Bianchi et al. (2007) confirm the pres- 
ence of the X-ray Baldwin effect in different samples of 
AGNs. Jimenez-Bailon et al. (2005) question the sig- 
nificance of the effect by finding only a 30% probability 
that the EW and the 2-10 keV luminosity are correlated. 
They and Jiang et al. (2006) both argue that the effect 
is mainly due to the radio-loud sources in the sample. 

The main problem is that the fitted correlation for 
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Fig. 8. — X-ray Baldwin effect: neutral Fe line equivalent width 
(EW) against the 2-10 keV luminosity. The solid line shows the 
best-fit to the AGN sample from Bianchi et al. (2007); the dotted 
lines represent the error on their fitted correlation. Filled circles 
show our RBS objects. 



the X-ray Baldwin effect comes mainly from AGNs with 
10 42 < L 2 -w/[ergs- 1 ] < 4 x 10 44 . Jimenez-Bailon et al. 
(2005) and Bianchi et al. (2007) note that more luminous 
radio-quiet QSOs with high S/N are needed for a solid 
analysis of the X-ray Baldwin effect. Our study adds 
three data points to this very poorly studied luminos- 
ity range, with EW= 70 - 140 cV ((EW) = 110). This 
is still consistent with the PG quasar average value of 
(EW) = 80 ± 20 cV with an dispersion of ct E w < 40 
cV (Jimenez-Bailon et al. 2005). Based on our 2-10 
kcV luminosities, we calculate the expected EW from 
the best- fit of Bianchi et al. (2007) (their Equation 1, 
sample of radio-quiet AGNs). The expected values of 
EW~ (31 — 40) ± 6 eV are significantly lower than what 
we measure in our RBS-QSOs (Figure 8). Except for 
RBS 1124 (XMM-Newton and Suzaku observations), the 
values disagree at the >90% confidence level. However, 
other studies also find much higher EWs (as predicted by 
the X-ray Baldwin effect) of -100 eV for L 2 -io > 10 4 ' 4 
erg s _1 (e.g., Chaudhary et al. 2010; Reeves & Turner 
2000; Page et al. 2004c). Figure 8 only contains objects 
from Jimenez-Bailon et al. (2005) that have a line energy 
consistent with E = 6.4 keV. Their sample and the sam- 
ple of Bianchi et al. (2007) have 25 PG quasars in com- 
mon. Although consistent within their uncertainties, the 
derived EW values for detected Fc lines by Bianchi et al. 
(2007) are systematically lower because, in contrast to 
Jimenez-Bailon et al. (2005), they fixed the line energy 
to E — 6.4 keV (S. Bianchi, private communication). 

Based on our study of RBS-QSOs, and on the Lx— EW 
relations of Chaudhary et al. (2010) and Jimenez-Bailon 
et al. (2005), we conclude that the X-ray Baldwin effect, 
the proposed anti-correlation between X-ray luminosity 
and EW of the narrow Fe line, is only present up to 
£2-10 — 10 44 erg s _1 . Above this luminosity thresh- 
old, there is convincing evidence for a flattening of the 
X-ray Baldwin effect or even a luminosity independent 
value of the equivalent width of (EW) — 100 eV. We 
hereby redefine the X-ray Baldwin effect not only as the 
proposed anti-correlation between Lx and EW but also 
as the flattening above L2-10 ~ 10 44 erg s _1 to cover 



a luminosity range as wide as possible. With current 
data samples, we are not able to tell if this flattening 
applies only to radio-quiet QSOs or for all QSOs. How- 
ever, the flattening would explain the discrepancy in the 
significance of the X-ray Baldwin effect when different 
luminosity ranges are explored, such as the very high lu- 
minosity range for PG quasars in Jimenez-Bailon et al. 
(2005).' 

The most fundamental question concerning the X-ray 
Baldwin effect is the physical cause why the behavior 
changes above L2-10 ~ 10 44 erg s _1 . Depending on the 
origin of the Fe line, different explanations are suggested. 
If the flattening or constant EWs above L2-10 ~ 10 44 
erg s _1 prove to be solid, either some physical process 
reaches its limit at this luminosity or two different mech- 
anisms/origins of the Fe line are causing the X-ray Bald- 
win effect. For our objects, the X-ray Baldwin effect is 
most likely not caused by the notion that an increase in 
X-ray luminosity progressively further ionizes the upper 
layers of the disk and therefore decreases the strength of 
the line (Ross & Fabian 2005) because the main contribu- 
tion to the observed Fe line does not seem to be located 
in the inner disk region. Other possible explanations are 
that with increasing luminosity, the inner radius of the 
torus surrounding the supermassive black hole increases 
up to a certain limit or the covering factor or thickness of 
the torus is decreasing (Page et al. 2004a). If the origin 
can be constrained to be in the torus, this luminosity- 
dependent covering factor may be directly linked to the 
decrease of the fraction of obscured AGNs (Ueda et al. 
2003). Both scenarios lead to a lower amount of repro- 
cessed radiation and therefore lower observed neutral re- 
flection and narrow Fe emission line strengths. However, 
the model still has to explain the high accretion rates 
found in RBS-QSO that seem to be sustained over a long 
period of time. 

10. CONCLUSIONS 

We select the 12 most luminous radio-quiet QSOs from 
the ROSAT Bright Survey (RBS) and present the first 
mini-survey of four of these sources. They each have 
£0.5-2 _ 10 45 erg s _1 , radio-to-optical flux densities 
R < 10, and allow us to study the accretion disk proper- 
ties at the very high end of the QSO luminosity function. 
These four objects have z = 0.21 — 0.55 and each have 
been observed with XMM-Newton; RBS 1124 has also 
been observed with Suzaku (Miniutti et al. 2010). Using 
emission lines from optical spectra, we estimate black 
hole masses of log(M BH /M Q )= 8.06 - 8.88 and Edding- 
ton ratios of iboi/^odd = 0.10 — 1.43 for these objects. 
Therefore, the outstanding luminosities are a combina- 
tion of high but not extreme black hole masses and Ed- 
dington ratios. The XMM-Newton 0.5-2 keV fluxes vary 
only ~5%-20% from the ROSAT fluxes, while the Suzaku 
and XMM-Newton observations of RBS 1124 show a 50% 
flux increase in only 13 months. The optical-to-X-ray 
spectral indices of our objects cover the same range found 
in other R OSA T selected AGN samples. 

All 0.2 12 keV XMM-Newton spectra show very ho- 
mogenous properties, implying that the physical proper- 
ties in the individual accretion disks are most likely very 
similar in this uniformly selected sample. The data of 
the three high S/N objects are well described by either 
(1) a model with two power laws or (2) an incident power 
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law and its reflection off the accretion disk very close to 
the black hole, but in all cases, we also have to model the 
presence of a narrow/barely resolved Fe line originating 
in neutral material light-days or further from the central 
X-ray source with an additional Gaussian line profile. 

Reflection models alone are not able to explain the 
strength of the Fe line. Two degenerate disk reflection 
model scenarios, which lie in seemingly very different 
areas of parameter space, are found for the high S/N 
spectra: (1) a disk with a very low ionization parame- 
ter (REFL1) with a reflection component which makes 
up only ~5%-20% of the total output or (2) a disk with 
a very high ionization parameter (REFL2-as almost ex- 
pected for very luminous QSOs) with a reflection com- 
ponent that accounts for ~35%-75% of the total output. 
Due to the weak contribution of the reflection compo- 
nent in REFL1, even the limited S/N RGS data do not 
allow us to favor either model. This degeneracy in the 
disk reflection model has to be broken by either more 
advanced models (maybe disk reflection models are even 
inappropriate for very luminosity radio-quiet QSOs) or 
extremely high-resolution spectra such as those which 
will be obtained by upcoming X-ray missions such as 
Astro-H (Takahashi et al 2008) and IXO (White et al. 
2009) to identify the physical properties of the accretion 
disk in the most luminous radio-quiet QSOs. Both sce- 
narios are consistent with the presence of fast rotating 
Kerr black holes. 

The soft excess in each source is well fitted by a power 
law, and its strength (2%-26% of the total 0.5-2 keV 
emission) is weaker than in lower luminous radio-quiet 
AGN/QSOs where the soft excess is typically 30% of the 
total 0.5-2 keV emission. No significant evidence for neu- 
tral absorption in excess of the Galactic absorption and 
for partial/full-covering absorption by ionized material is 
found. The RGS data confirm the extreme smoothness 
of the soft excess and the absence of ionized absorbers. 

The major fraction of the Fe line emission in our three 
high S/N RBS-QSOs has to come from lowly-ionized or 
neutral material located only a few light days (or fur- 
ther) out, but a superimposed very weak relativistically 
broadened Fe line from the inner disk cannot be ruled 
out (90% confidence level upper limits from fitting laor 
components: ~0. 14-0.40 keV). It remains unclear how 
a considerable fraction of the line-emitting material can 
avoid becoming highly ionized so close to the central X- 
ray source of the most luminous QSOs. 
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jects, we can confirm the stacking technique assumption 
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interpretation of the low S/N spectra of QSOs with sim- 
ilar luminosities at higher redshifts, which are routinely 
found by XMM-Newton and Chandra surveys. Knowl- 
edge about intrinsic spectral components and the scatter 
in their parameters values is important e.g., when esti- 
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the best-fit results for a model with two power laws plus 
a Gaussian Fe line profile in Section 8). 
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